Comparison among Various Expressions of Complex Admittance 
for Quantum System in Contact with Heat Reservoir 



Mizuhiko Saeki 1 ' 3 , Chikako Uchiyama 2 , Takashi Mori 1,3 and Seiji Miyashita 1,3 

1 Department of Physics, The University of Tokyo, 
7-3-1, Hongo, Bunkyo-ku, Tokyo 113-0033, Japan 

2 Faculty of Engineering, University of Yarnanashi, 
4-3-11, Takeda, Kofu, Yarnanashi 400-8511, Japan 

3 CREST, JST, 4-1-8, Honcho, Kawaguchi, Saitama 332-0012, Japan 

(Dated:) 

Relation among various expressions of the complex admittance for quantum systems in contact 
with heat reservoir is studied. Exact expressions of the complex admittance are derived in various 
types of formulations of equations of motion under contact with heat reservoir. Namely, the complex 
admittance is studied in the relaxation method and the external-field method. In the former method, 
the admittance is calculated using the Kubo formula for quantum systems in contact with heat 
reservoir in no external driving fields, while in the latter method the admittance is directly calculated 
from equations of motion with external driving terms. In each method, two types of equation of 
motions are considered, i.e., the time-convolution (TC) equation and time-convolutionless (TCL) 
equation. That is, the full of the four cases are studied. It is turned out that the expression of 
the complex admittance obtained by using the relaxation method with the TC equation exactly 
coincides with that obtained by using the external-field method with the TC equation, while other 
two methods give different forms. It is also explicitly demonstrated that all the expressions of the 
complex admittance coincide with each other in the lowest Born approximation for the system- 
reservoir interaction. The formulae necessary for the higher order expansions in powers of the 
system-reservoir interaction are derived, and also the expressions of the admittance in the n-th 
order approximation are given. By transforming inverse-temperature-integrals into time-integrals, 
the admittances are also given in the formulae with time-integrals alone. To characterize the TC 
and TCL methods, we study the expressions of the admittances of two exactly solvable models. 
Each exact form of admittance is compared with the results of the two methods in the lowest Born 
approximation. It is found that depending on the model, either of TC and TCL would be the better 
method. 



PACS numbers: 

I. INTRODUCTION 

The complex admittance is a fundamental quantity to 
describe response of systems to an external field. Since 
the magnetic resonance has been discovered [1], exper- 
imental and theoretical studies have been developed si- 
multaneously [2-4]. For example, the magnetic resonance 
(MR) has been one of the most active fields to study the 
quantum response. Kubo and Tomita gave a unified view 
point to express the response function by the autocorre- 
lation function [5], which was extended to the general 
formula of linear response [6] . 

Although the Kubo formula [6] gives a general form of 
the complex admittance, there are several types of formu- 
lation to take into account effects of the interaction of a 
quantum system with its heat reservoir. When we study 
the time dependence of a physical quantity A under a 
time-dependent external field F conjugate to a physical 
quantity B, the response function can be expressed in 
terms of a time correlation function of quantities A and 
B, i.e., (A(t)B). In the simplest treatment, we give the 
time evolution by the pure quantum dynamics given by 
Hamiltonian Hs of the system. In this case, the complex 
admittance xab (w) is given by an ensemble of delta func- 



energy levels E^ and Ej of the system. 

In many-body systems, the energy levels form a con- 
tinuous spectrum, and also the complex admittance has 
continuous forms. There, line shape of the spectrum (the 
imaginary part of the complex susceptibility) has an in- 
trinsic line width which is attributed to the energy struc- 
ture of the system. In such systems, effects of interac- 
tions on the spectrum line shape have been developed 
[7], and various properties have been clarified. For ex- 
ample, Nagata and Tazuke found the shift of resonance 
peak due to dipole-dipole interaction as a function of rel- 
ative angle between the external field and the lattice axis 
of one-dimensional magnets [8, 9]. Recently, studies to 
evaluate the susceptibility from microscopic Hamiltonian 
have been proposed, and new aspects of the resonance 
have been discussed [10, 11]. 

On the other hand, the width of the line shape also 
appears due to the interaction of system with its heat 
reservoir. Because of the developments of experimental 
methods, real time quantum dynamics has become ac- 
cessible in microscopic quantum systems, e.g., the mag- 
netization dynamics of single molecule magnets [12-20] 
and also microscopic circuits which manipulate the de- 
gree of freedom of qubit [21]. In these systems, it has 
been pointed out that the dissipation plays important 
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been discussed in detail [24]. It is interesting problems 
to study effects of the dissipation on the complex admit- 
tance for those systems. 

If we consider the total system which consists of a sys- 
tem of interest and of its heat reservoir, the formulation 
of a pure quantum dynamics works in principle. How- 
ever, it is difficult in practice, and we introduce a kind 
of dissipative dynamics instead of the quantum dynamics 
for the total system. Usually we introduce the dissipa- 
tive dynamics by coarse-graining of the heat reservoir. 
The dissipative dynamics is expressed by the so-called 
quantum master equation [25-32]. 

In order to take into account the dissipative dynam- 
ics, we need to obtain the time evolution function for 
dissipative systems. The phenomenological treatment 
of the dissipation effects had been done by using the 
Bloch equation including disipation terms with the re- 
laxation times T\ and T 2 [3] . The microscopic treatment 
was proposed using the projection operator techniques 
by Nakajima [25], Zwanzig [26] and Mori [33]. They 
obtained the equations of motion in the non-Markovian 
or "time-convolution (TC)" integral forms with memory. 
However, in order to treat the memory term in a com- 
pact way, the so-called "time-convolutionless (TCL)" for- 
malism! has been introduced [28, 34, 35]. Recently, we 
have studied a complete expression of the complex admit- 
tance in the dissipative environment in the TC formalism 
[36, 37]. There, we have used a standard equation of the 
reduced density operator, but we have pointed out that 
the initial correlation between the system and the heat 
reservoir can not be ignored in general when we treat the 
non-Markovian case. 

Moreover, one of the authors has proposed the formal- 
ism called the "TCLE method" in which the complex ad- 
mittance is directly calculated from time-convolutionless 
(TCL) equations of motion with external driving terms 
[38-44], by generalizing the method of Argyres and Kel- 
lcy in which the admittance is calculated from time- 
convolution (TC) equations of motion with external driv- 
ing terms in the Markovian approximation [45]. Here- 
after, we call the method in which the complex admit- 
tance is directly calculated from equations of motion 
with external driving terms, the "external-field method" , 
while we call the method in which the Kubo formula is 
calculated for systems with no external driving fields, the 
"relaxation method" . Thus, we have the four types of for- 
malism, that is, the relaxation method or the external- 
field method with TC or TCL equation of motion. The 
external-field method is considered to give the same re- 
sults as those obtained using the relaxation method in 
principle. They must also give the same results because 
the original equation is the same. However, the treatment 
of the contact with the heat reservoir can be simplified 
in the external-field method. In practical calculations, 
we need to approximate perturbatively in order of the 
strength of the interaction with the heat reservoir. If we 
truncate the perturbation, the higher order corrections 

i_ _ _i*£r j_ • _ r r !• 



In the present paper, we consider a system interact- 
ing with a heat reservoir in an external driving field, and 
derive forms of the complex admittance using the four 
types of formalism. Then, we investigate relations among 
them. Namely, we examine the forms of admittance ob- 
tained by using the following four methods. 

(1) The relaxation TC method (the relaxation method 
with the TC equation) 

(2) The relaxation TCL method (the relaxation method 
with the TCL equation) 

(3) The TCE method (the external-field method with the 
TC equation) 

(4) The TCLE method (the external-field method with 
the TCL equation) 

We also derive the formulae necessary for the higher or- 
der expansions in powers of the system-reservoir interac- 
tion, and give the expressions of the admittance in the 
n-th order approximation. Moreover, in order to dis- 
cuss the truncation effect in the four methods, we study 
the two exactly solvable models, and compare each exact 
form of admittance with the results obtained by using 
the above methods in the lowest Born approximation for 
the system-reservoir interaction. 

In Section II, we survey the derivation of the Kubo 
formula, and explain the basic difference between the re- 
laxation method and external-field method. In Section 
III, we derive forms of the complex admittance obtained 
by using the above four methods. In Section IV, we inves- 
tigate relations among forms of the admittance obtained 
by using the four methods in the lowest Born approxi- 
mation for the system-reservoir interaction. In Section 
V, we derive the formulae necessary for the higher order 
expansions in powers of the system-reservoir interaction. 
In Section VI, we study the two exactly solvable models 
to discuss the truncation effects in different methods. In 
Section VII, we give a short summary and some conclud- 
ing remarks. 



II. MODEL 

We consider a quantum system interacting with a 
quantum heat reservoir in an external driving field. We 
take the Hamiltonian TLj{t) of the total system as 

W-r(i) = T~ts + %r + TtsR + T~ted(t) 

= Ho + H SR + H ed (t) = n + H ei (t), (1) 

where Tis is the Hamiltonian of the quantum system, 7i R 
is the Hamiltonian of the heat reservoir, Hsr is the inter- 
action Hamiltonian between the system and heat reser- 
voir, and 7i e d(t) is the interaction Hamiltonian of the 
quantum system with the external driving field. Here, we 
express the Hamiltonian of the system and heat reservoir 
with interaction by TL = Hs + Hr + Hsr, and the unper- 
turbed Hamiltonian, i.e., the Hamiltonian of the system 
and heat reservoir without interaction by Ho = Hs + Hr. 
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the Liouville equation 

(d/dt) p T {t) = -(i/h) [H T (t) , p T {t) ] = -i L T (t) pr(t), 

(2) 

which decides the dynamics corresponding to the Hamil- 
tonian Hi(t), where the Liouvillian Ljit) corresponds to 
the Hamiltonian Hj(t). Hereafter, a Liouvillian corre- 
sponding to a Hamiltonian, say H, is denoted as L, i.e., 
LA=[H,A]/h. 

We assume that the system and heat reservoir are in 
the thermal equilibrium state Pte at temperature T ini- 
tially, i.e., before the external driving field is turned on, 
where Pte is the thermal equilibrium density operator for 
the system and reservoir and is given by 

p TE = cxp(-/3H)/Trexp(-/3H), (3) 

with (}=\/{k%T). Here, notation Tr denotes the trace 
operations in the system and reservoir spaces, i.e., 
Tr = trtr R , where notations tr and tr R denote the trace 
operations in the system and reservoir spaces, respec- 
tively. If we consider only the system, the density oper- 
ator for the system alone is given by 

p s = exp(-/3W s )/trexp(-/3W s ). (4) 

The density operator for the heat reservoir alone is given 

by 

(OR = exp(-/3W R )/tr R exp(-/3W R ). (5) 

It should be noted that the equlibrium state of the system 
is given by 

Pa = tr R ^ = tr R exp(- /? H)/ Tr exp(- H) , (6) 

which includes the interaction between the system and 
reservoir and is different from p s . 

For the convenience of expressions, we define the nota- 
tion (• • • ) R = tr R (- • • p R ), and renormalize the Hamiltonian 
Hs of the system and the system-reservoir interaction 
Hsr as 

Hs — > H s + (W SR ) R =► H s , (7a) 
H S r — > TC S r - (Hsk)r =>■ H S r- (7b) 

Hereafter, we use Hs and Hsr for the Hamiltonians renor- 
malized by (7), and then, we have (7^sr)r = 0. 

We now take the interaction of the system with the 
external driving field as [6, 45] 

H ed (t) = -$>^.(i) = -Y^ A i F Me- iu \ 

3 3 w 

= ^HedHe- ,ui , (8) 

where Aj is the physical quantity conjugate to the force 
Fj(t) which is a c-number function of the time t, and the 

• . i . r___ - _ . 



A. The Kubo formula 

We first survey the derivation of the Kubo formula 
briefly [6]. The time evolution for the density operator 
of the total system is given by Eq. (2). The external 
driving field is assumed to be turned on adiabatically at 
the initial time t = to which is infinite past (to = — oo) 
in the Kubo theory [6]. Expanding the density operator 
Pi(t) in powers of the external driving field as 

pr(t) = pro(t) + pn{t) + pn(t) + • ■ ■ , (9) 

the zeroth-order part pio(i) and first-order part pn(t) 
satisfy the following equations and initial conditions : 

^ Pm(t) = - iLp w {t) ; Pm(to) = Pie , (10) 

^ Pti (t) = - iL pn (t) - iL ed (t) p T0 (*) ; fhi (to ) = . 

(11) 

Equation (10) has, by virtue of the thermal equilibrium 
density operator (3), the solution 

p T0 (t) = exp{-iL (t - t )}pm(t ) = Pte , (12) 

by which Eq. (11) can be formally solved as 

(>n(t) = -i dr exp{- i L ■ (t - t)} L ed (r) pit . (13) 

Then, the first-order part in powers of the external driv- 
ing field for the expectation value of a physical quantity 
Ai of the system, can be described as 

Tr Ai pn{t) 

= \ E E f dT Tr A * e ~ l L (t ~ T) . m *~ ' " T > 

= ^EE f ^ dTTvAie-^^A^^FM 

j w 

xe tu ' T - tu ' t . (14) 

The complex admittance Xij(uj) is defined in the limit 

to — > — oo, as [6] 

Tr A l pn = Xii M Fj M e~ iut , (15) 

3 w 

(t — > - oo) 

and takes the expressions [6] 
X«(w) = % -j™dtTrA i e- iLt [A j , p TE ] e ^*- e *, (16a) 




dtTrAi [A^(—t), Pte] e iujt ~ et , (16b) 
dtTtA^(t)[A h pTE}e iu,t - et , (16c) 



= -/ dtTviA^A^p^e^-^, (16d) 
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which are called the Kubo formula, where Here, 
A^(t) is the Heiscnbcrg operator defined by 

A w (t) = exp(iLt)A = exp(iHt/ K)Aexp(-iHt/ h). 

(17) 

In the definition (15) of the complex admittance, the 
summation for lo is over the same frequencies as 
those in (8). In the Kubo formula, A?(t) and A?(t) denote 
the time evolutions of the physical quantities A^ and Aj . 
Suzuki and Kubo discussed the linear response by deriv- 
ing the time evolution of the physical quantity A t using 
the stochastic dynamics [46] . In the Kubo theory [6] , the 
complex admittance is given by the time correlation func- 
tion for the system with no external driving fields. We 
call the method that evaluates the complex admittance 
for a quantum system in contact with its heat reservoir 
using the Kubo formula, the "relaxation method" . 

We next mention the basic difference between the re- 
laxation method and external-field method in the follow- 
ing subsections. 



B. The relaxation method 

The relaxation method is the method in which com- 
plex admittances are calculated using the Kubo formula 
for a quantum system in contact with its heat reservoir. 
For a quantum system in contact with its heat reservoir, 
the operator with a physical quantity A of the quantum 
system : 

A(t) = e- iLt [A,p TE ] =e- iHt l h [A,f^]e iUt l h , (18) 

satisfies the equation of motion 

(d/dt) A(t) = -(i/h)[H, A(t) } = -iL A(t). (19) 

In order to eliminate irrelevant variables associated with 
the heat reservoir, we introduce the time-independent 
projection operators V and Q, which are taken as 



V = Pr tr R , 



Q = \-V 



tr R p R = l, (20) 



where p R is the equilibrium density operator (5) of the 
heat reservoir alone. Applying the projection operators 
V and Q to Eq. (19), we have the coupled equations : 

(d/dt)VA(t) = -iVLVA(t)-iVLQA(t), (21a) 
{d/dt) QA(t) = - i QLV A(t) - i QLQA(t). (21b) 

In the relaxation method, complex admittances are cal- 
culated by deriving the equation of motion for tr K A(t) 
from the above two equations (21) and by substituting 
its solution into the Kubo formula (16). 



equation of motion including external driving terms for 
the total system, which includes the external driving 
field. The density operator pj(t) of the total system sat- 
isfies the Liouville equation (2). We assume that the 
external driving field is turned on adiabatically at the 
initial time to = 0, though the initial time is infinite past 
(to = — oo) in the Kubo theory [6]. Expanding the den- 
sity operator pi(t) in powers of the external driving field 
as in (9), the zeroth-order part pm(t) and first-order part 
Pn{t) satisfy the following equations and initial condi- 
tions similar to (10) and (11) : 

^ Pm(t) = - iL p m (t) ; /Oto(0) = Pte, (22) 

^PTi(t) = -iLf>n(t) - iL ed (t)p m (t) ; pri(0) = 0. 

(23) 

Applying the projection operators V and Q defined by 
(20) to Eq. (23), we have the coupled equations : 

{d/dt)V Pri(t) = - iVLVpnit) - iTLQp-nit) 

-iVL ed (t)p TE , (24a) 

(d/dt) Qpn(t) = - iQLVpn{t)- iQLQp T1 {t) 

-iQL ed {t) PlE , (24b) 

where we have used the solution of Eq. (22) : 

Pro{t) = cxp( i Lt) Pto(0) = Pte ■ (25) 

In the external-field method, complex admittances are 
calculated by deriving the equation of motion for 
tr R pii (t) from the above two equations (24) and by solv- 
ing it in the limit i^oo, which corresponds to the as- 
sumption that the initial time is infinite past (to = — oo) 
in the Kubo theory [6] . 

In the external-field method, it should be noticed that 
the time evolution of the external driving force is included 
in the projection procedure. Thus, the projection proce- 
dure is performed for the equation (23) of motion includ- 
ing external driving terms in the external-field method, 
while the projection procedure is performed for the equa- 
tion (19) of motion including no external driving terms 
in relaxation method. 



III. COMPLEX ADMITTANCE 

In this section, we derive forms of complex admit- 
tance Xij( w ) f° r t ne tw° methods mentioned in the pre- 
vious section, i.e., the relaxation method and external- 
field method, and for the two types of equation of mo- 
tion, i.e., the time-convolution (TC) equation and time- 
convolutionless (TCL) equation. 



The external-field method 



The Relaxation TC method 



The external-field method is the method in which com- 



We first derive the form of complex admittance by the 



calculated by solving the "time-convolution" (TC) equa- tion of motion for systems with no external driving fields. 

I 

Equation (21b) has the formal solution 

QA(t) = -i[ rfr exp{- iQLQ(t- r)} QLV Mr) + exp(- iQLQt) QA(0), (26) 
Jo 

which has a form of the time-convolution (TC). By substituting (26) into (21a), the "time-convolution" (TC) equation 
of motion for VA(t) can be obtained as 

(d/dt)VA(t) = -iVLVA(t) - [ dr VL exp{- iQLQ(t-r)}QL VA(t) - i VL exp(- iQLQt)Q A(0). (27) 

Jo 

By virtue of the projection operators (20), the above equation is reduced to the TC equation of motion 

(d/dt) a(t) = -iL s a(t) + C(t , {a}) + I(t), (28) 
for the reduced operator a(t) defined by 

5(t) = tr R A(t) = tr R exp(— iLt)[A, p TE ] , (29) 
where the collision term C(t, {a}) and the inhomogeneous term I{t) are, respectively, given by 

C(t,{d}) = - [ dr tr R _L SR exp{— i Q L Q (t — t)} L sr p R 5(t), (30) 
Jo 

I(t) = -«-tr R L S RCxp(-zQLQf)Q[A, p TE ]. (31) 

The inhomogeneous term I{t) represents the effects of initial correlation of the system and reservoir, because if the 
initial state pte given by (3) is approximated to the decoupled one poPr, i-e., Pte = Po Pr, the term I(t) vanishes since 
Q [A, pte] = Q [A, po] p R = 0. Performing the Fourier-Laplace transformation for Eq. (28), we have 

-o(0) -iwo[w] = -iL s a[uj] + C[w] a[w] + I[u], (32) 

where d[u], C[u>] and I[oj] are given by 

d[iv] = / dt d(t) exp(i uo t — e t) = I dt tr R exp(— i L t) [A , p-^] exp(i u)t — et)\ _^ . , (33) 
Jo Jo e ^ 

/•OO 

C[w] = - dt(L sn exp{-iQLQt)L sn ) R exp{iujt-et)\ e ^ +0 , (34) 

J 

I[ui\ = -iJ efttr R L SR exp(-i QL Qt) Q [A, p TE ] cxp(iujt - e£)| e ^ +0 • (35) 
Considering that 

a(0) =tr R i(0) =tr R [A, p TE ] = [A,p ], (36) 
Eq. (32) can be formally solved as 

5[w] - (i (Ls - w) - CM) _1 { [ A , po ] + I[w] } . (37) 
Substituting (37) into the Kubo formula (16), the admittance takes the form 

X^j C (uj) — — / dt Tr Ai exp(— i Lt) [Aj , p TE ] exp(i tot) = — / dt tr A t dj(t) exp(iwi) = — tr Ai dj[oj], 
n J n Jq h 

= trA * 777 rTTVTT^r ll > ^ + ^M}, (38) 
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which is a general form of admittance derived using the relaxation TC method under thermal equilibrium initial 
conditions, where Dj [u] comes from the inhomogeneous term I(t) in Eq. (28), represents the effects of initial correlation 
of the system and reservoir and is given by 

i 1 f°° 

Dj [u] = - Ij [u] = - J dt tr R L SR exp(- i Q L Q t) Q [A 3 , p TE ] exp(i u t - e i) \ e ^ +Q . (39) 

Expanding Eq. (28) up to second order in powers of the system-reservoir interaction Hs R , it reduces to 

(d/dt)a(t) =-iL s a{t)+CW(t, {d})+I (2) (t), (40) 
where the collision term C^ 2 \t, {a}) and the inhomogeneous term I^ 2 \t) are, respectively, given by 

(t, {5}) = - / dr (L SR cxp{- i L (t - r)} L SR > o(r), (41) 
Jo 

[ & 

(t)=i dp' tr R L SR exp(- i L t) [A , p s p R W SR (- i h (?) ] . (42) 
Jo 

Here, Hs R {t) is defined by 

7i SR (t) = cxp(z L t) H SR = cxp(i H t/K) U SR exp(- i H t/h). (43) 

Performing the Fourier-Laplace transformation for Eq. (40) and substituting its formal solution into the Kubo formula 
(16), the admittance takes the form 

- ^ i{La _^_ CW[u] {i [Aj , Ps + ] + D?[u]}, (44) 

which is a general form of admittance derived using the relaxation TC method in the lowest Born approximation for 
the system-reservoir interaction, where pg 2 \ C l ^ 2 - ) [a;] and L'j 2 ''^] are, respectively, given by 

p ( o 2) =Ps d/3 2 {tv R p R H SR (-ihp 1 )H SR (-ihp 2 )-Trp s p R H SR (-ihp 1 )H SR (-ihp 2 )}, (45) 

J ^0 

(-00 rOO 

C (2) [uj} = - rfr(L SR cxp(-2 J LoT) J L SR ) R cxp(zwT) = - / dr (L SR L sr (- t)) r cxp{i (u - L s ) t}, (46) 
Jo Jo 

bf [u] = - lf ] [w] = — J dr J d(3' tr R L SR exp(- zL t)[ A 3 , p s p R H SR (- ih[3')] exp(i w r). (47) 



The term [u>] represents the effects of initial correlation of the system and reservoir. The above result (44) can 
be obtained by expanding Eq. (32) up to second order in powers of the system-reservoir interaction too. By using 
this formula, we can make a unified numerical method to treat systems with large number of spins, and can study 
relations between the spin-spin interaction and the noise effects on the line shape of magnetic resonace [37] . 

B. The Relaxation TCL method 

One of the authors has derived the form of complex admittance by the "relaxation TCL method" in which the Kubo 
formula is calculated by solving the "time-convolutionless" (TCL) equation of motion for systems with no external 
driving fields [43, 47]. We here survey the derivation of the form of complex admittance by the relaxation TCL method 
for comparison. In order to renormalize the time-convolution (TC) in the formal solution (26) of Eq. (21b), we write 
the formal solution of Eq. (19) in the form 

A(t) =exp{-iL-(t-T)}A(r) or A(t) = cxp{i L ■ (t - r)}A(t), (48) 

which is substituted into (26) to give 
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where 9(t) is defined by 

9{t) = { l + l J Q rfrexp(-iQLQr)QL SR Pexp(iir)} \ (50) 

By substituting the expression (49) for QA(t) into (21a), the "time-convolutionless" (TCL) equation of motion for 
VA(t) can be obtained as 

[d/dt)VA{i) = -iVLVA{i) - iVL{9(t) - l}VA(t) - iVL9(t) exp(- iQLQt) QA(0). (51) 

By virtue of the projection operators (20), the above equation is reduced to the TCL equation of motion for the 
reduced operator a(t) as 

(d/dt) a(t) = -iL s a(t) + C{t) a(t) + I(t), (52) 
where the collision operator C(t) and the inhomogeneous term I(t) are, respectively, given by 

C(t) = - i ■ tr R L SK {6{t) - 1} PK ee - i { L SR {0(t) - 1} } R , (53) 
I(t) = —i • tr R L SR #(i) exp(— i QL Qt) Q[A, p TE ] . (54) 

Equation (52) corresponds to the TCL equation of motion for the reduced density operator of the system [28]. The 
inhomogeneous term I(t) represents the effects of initial correlation of the system and reservoir, because if the initial 
state Pte [(3)] is approximated to the decoupled one poPr, the term I(t) vanishes. Equation (52) has the formal 



(55) 



solution 

a(t) = cxp^{-iL s t + drC(r)} a(0) + drcxp^j- i L s ■ (t - t) + J dsC(s)} I(t). 

Substituting (55) into the Kubo formula (16), the admittance takes the form [43] 

if 00 if 00 
X*] CL (w) = — / dt Tr Ai exp(- i L t) [ A i , p TE ] exp(z ut) = - dttrAi hj (t) cxp(i oj t), 
ft Jo ft Jo 

= -J dt tr Ai cxp^|i(w - L s )t + J drC(T)j [Aj , p ] 

+ ~ftl dt I rfTtr ^ exp ^{*( w_jLs )( <_T ) + / dsC(s)} Ij(T)exp(iujT), (56) 

which is a general form of admittance derived using the relaxation TCL method under thermal equilibrium initial 
conditions, where Ij(t) is equal to I(t) with Aj in place of A. The first term of the above admittance (56) is the 
admittance obtained using the conventional relaxation method under the decoupled initial condition /5t(^o) = Pte =^ 
Po Pr, and its second term comes from the inhomogeneous term I{t) in Eq. (52) and represents the effects of initial 
correlation of the system and reservoir. 

Expanding Eq. (52) up to second order in powers of the system-reservoir interaction Hsr, it reduces to 

(d/dt) a(t) = -iL s 5(t) + C {2) (t) d(t) + 7 (2) (t), (57) 

with 

C (2) (*) = - / rfr(L SR exp(-iioT)L SR exp(ii T)) R = - / dr (L S rL S r(- r)> , (58) 
Jo Jo 

I M(t) = lV)(t)=i d/3'tr R L SR exp(-ii t) [A, p s PrH sr {- i h p')] . (59) 



Solving Eq. (57) formally and substituting its solution into the Kubo formula (16), the admittance takes the form 
[42, 43] 

X*f L H = l fri°° dt tr Ai Cxp - {' ("-Ls)t + f*dT (r) } [Aj , p s + p[ 2) ] 

+ % ft\ ^ I dTtrj4 ' ex P-{ l ( w_is )(^ T )+ / dsC (2) (s)|/j 2) (T)cxp(iwT), (60) 

which is a general form of the admittance derived using the relaxation TCL method in the lowest Born approximation 
for the system-reservoir interaction, where lf\t) is equal to I^(t) with Aj in place of A. The second term of 
the above admittance comes from the inhomogeneous term I^ 2 \t) in Eq. (57) and represents the effects of initial 
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C. The TCE method 



We next derive the form of complex admittance by the "TCE method" in which complex admittances are directly 
calculated by solving the time-convolution (TC) equation of motion including external driving terms. Equation (24b) 
has the formal solution 

Qp T1 (t) = -i[ dr exp{— i QL Q(t— t)} \QLV p^\(r) + Q L ed (r) Pte}, (61) 
Jo 

which has a form of the time-convolution (TC). By substituting the time-convolution expression (61) for Qpnif) into 
(24a) directly, the time-convolution (TC) equation of motion for VpTi(t) can be obtained as 

(d/dt) Vpji(t) = — i VL V pj\ (t) — VL [ dr exp{- i Q L Q(t - t)} QLVp^r) 

Jo 

- iVL ed (t)p TE -VL dr cxp{— i Q L Q (t — t)} Q L ed (r ) p TE . (62) 
Jo 

By virtue of the projection operators (20), Eq. (62) is reduced to the TC equation of motion for p\(t) [=tr R( o T1 (t)] 
[45]: 

(d/dt) Pl (t) = - i L s Pl (t) + C(t, {pi})-i L ed (t) p Q + D(t). (63) 
On the right-hand side of Eq. (63), the collision term C(t, {pi}) and the inhomogeneous term D(t) are given by 



c ( t APi}) = - / dr tr R L SR exp{- i QL Q (t - r)} XsrPrPi (t) = - / dr ( L SR exp(- i QL Qt) L SR ) R p 1 (t - r), 
Jo Jo 

(64) 

D(t) = - I dr tr R L SR cxp{- i Q L Q (t - r)} L ed (r) Q p TE = - f dr tr R L SR exp(- i QL Qt) L ed {t - r) Qpr E . (65) 
Jo Jo 



The inhomogeneous term D(t) is called the "interference term" in the TC equation, and represents the effects of initial 
correlation of the system and reservoir, because if the initial state p TE [(3)] is approximated to the decoupled one PoPr, 
the interference term D(t) vanishes. We now take into account the interaction of the system with the external driving 
field, which is turned on adiabatically at the initial time t = 0, to be given by the Hamiltonian (8). The admittance 
Xij(oj) for the physical quantities Ai and Aj of the system is defined in the limit t — > oo, as 

T±A i p T1 (t)=trA i p 1 (t) = J2J2x i i^ F ^ e ~ i " t > (t-oo) (66) 

j u> 

which corresponds to (15). We assume that the correlation times of the heat reservoir are finite, that the integrands 
of (64) and (65) vanish for finite r in the limit t^oo, and that Eq. (63) has a stationary solution of the form 



Pi 



(t)=J2piHe-^ t . (t-oo) (67) 



The summations f° r u m (66) and (67) are over the same frequencies as those in (8). Then, P\[lo] satisfies the 
equation 

- i lj pi[u>] = - i LsPiM + C[u)} pi[uj} - i L ed [uj] p + D[lu] , (68) 
which has the formal solution 

pM - (liLs-^-Cl^y^-iL^po + DlLu]}, (69) 
where C[w] is the collision operator given by (34), and D[u] is the interference term defined through the relation 

D(t) = y~] D[lj] e~ iuJt , (t^oo) (70) 
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and takes the expression 

poo 

D[uj] = - dr tr R L SR cxp(- i QL Qr) L ed [uj] QpT E exp(jwT). (71) 
Jo 

In the relation (70), the summation ^ u for ui is over the same frequencies as those in (8). By substituting the formal 
solution (69) into (67) and by using the definition (66) of the admittance, we can obtain 

*«' M - " A - TurdpcR is w • ""i + ^4 < 72 » 

which is a general form of the admittance obtained using the TCE method, where Dj [uj] is the interference term given 
by (39) and represents the effects of initial correlation of the system and reservoir. We notice that the admittance 
(72) obtained using the TCE method has exactly the same form as the admittance (38) obtained using the relaxation 
TC method. 

Expanding Eq. (68) up to second order in powers of the system-reservoir interaction 7i SR , it reduces to 

- i uj Pl [uj] = - i L s Pl [to] + C ( V [ u ] Pl [ u ] - i L ed [uj] ( Ps + Po 2) )+D& [uj] , (73) 
where the collision operator [uj] is given by (46) and the interference term [to] is given by 

/■oo r[3 

DW[w}= dr dl3'tr K L SR exp{-iLoT)L ed [w} PsPK H S n(-ihf3')eMiuT), (74) 
Jo Jo 

which is shown in Appendix A to coincide with 

/■OO PT 

D^[w]=i dr dstr R I/ S Rexp(- iL s) L ed [oj] exp{i Lq (s - t)} L S r P s p R exp(i u>s). (75) 
Jo Jo 

Substituting the formal solution of Eq. (73) : 

Pl [ U ] = (i (L s -uj)- [uj] ) - 1 { - i L ed [cu] ( Ps + p[ 2) ) + [uj] } (76) 
into (67) and using the definition (66) of the admittance, the admittance takes the form 

Xlf M - trJU i{La _^_ cW[u] {i [A 3 , Ps + Po 2) ] + ^M}, (77) 

which has the same form as the admittance (44) obtained using the relaxation TC method in the lowest Born 
approximation for the system-reservoir interaction. Here, the interference term D - [co] is given by (47), is shown in 
Appendix A to coincide with 

Df ) [uj] = --J dr J dsti K L SR exp(-iL s)[A j ,exp{iL (s-T)}L SRPspR ]exp(iujs), (78) 

and represents the effects of initial correlation of the system and reservoir. The above result (77) can be obtained by 
expanding the TC equation (63) up to second order in powers of the system-reservoir interaction too. In the TCE 
method, the admittance (77) is obtained by the lowest Born approximation of the equation (63) or (68) including 
external driving terms, while in the relaxation TC method, the admittance (44) is obtained by the lowest Born 
approximation of the equation (28) or (32) including no external driving terms. It should be noticed that although 
the equations in which the perturbation is truncated arc different from each other in the two methods, the expressions 
of the obtained admittances coincides with each other. 



D. The TCLE method 



One of the authors has proposed and has studied the "TCLE method" in which the complex admittance is directly 
calculated by solving the time-convolutionless (TCL) equations of motion including external driving terms [38-44]. 
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derive the exact form of complex admittance by the TCLE method and survey that method. In order to renormalize 
the time-convolution (TC) in (61), we write the formal solution of Eq. (23) in the form 

Pti(t) = exp{iL(t-<r)}p T i(i) +i J ds exp{i L (s - r)} L ed (s) p 1E , (79) 
where we have used the solution (25) of Eq. (22). Substituting (79) into (61) and solving it for Qpt\{t), we have 

QpTi(t) = {6(t)-l}Vp T1 (t)-i6(t) f dTexp{- i QL Q(t - t)} QO^ 1 ^) L ed (r) p TE , (80) 

Jo 

where 6(t) is given by (50). By substituting the above expression for Qpri(t) into (24a), the time-convolutionlcss 
(TCL) equation of motion for Vpri{t) can be obtained as 

(d/dt)Vp T1 (t) = - iVLVp T1 (t) -iVL{6{t) - l}V Pn {t) -tVL ed (t)p TE 

- [ dTVL0(t)exp{-iQLQ(t-T)}Q0- 1 (T)L«(T) f > ni . (81) 
Jo 

By virtue of the projection operators (20), Eq. (81) is reduced to the TCL equation of motion for pi(t) [= tr R pj\ (t)] 
[38, 39, 43] : 

(d/dt) Pl (*) = -* L s Pl (t) + C(t) Pl (t) - i L ed (t) Po + D(t). (82) 

On the right-hand side of Eq. (82), C(t) is the collision operator given by (53), and D{t) describes the effects of 
interference between the external driving field and heat reservoir, is called the "interference term" and take the form 

D(t) = - J dr tr R i S R 1 _ exp(- iQLQr) L ed (t - r) Q p 1E 

+ J dr tr R L SR exp(- i QL Qr)E(t - r) L ed (t - r) p TE , (83) 

with £(t) defined by 

Z(t) = l-Q- 1 (t) = -i f dr exp(— i QL Qt) Q L SR V exp(i L t). (84) 
Jo 

The first term of D(t) given by (83) is called the "first interference term" and represents the effects of initial correlation 
of the system and reservoir, because if the initial state p? E [(3)] is approximated to the decoupled one PoPr, this term 
vanishes. But, the first term of D(t) [(83)] includes the memory effects and is different from D{t) given by (65). The 
second term of D(t) is called the "second interference term" and represents the memory effects, which are the effects 
of collision of the system with the heat reservoir, because this term is derived by renormalizing the time-convolution 
in (61). We assume that the correlation times of the heat reservoir are finite and that Eq. (82) has a stationary 
solution of the form (67). Then, pi[u>] satisfies the equation 

- i wpi [u>] = -iL s p 1 [lu] + C pi [to] - i L ed [uj] p + D[to] . (85) 

On the right-hand side of Eq. (85), C is the collision operator given by 

C = C(oo) = - i ■ tr R L SR E (1 - S)- Vr = - * ■ <£sr £ (1 - S)" 1 ) R , (86) 

and D[to] is the interference term defined through the relation 

D(()^J}He~ iut , (t-foo) (87) 

and is shown in Appendix B to take the expression 
,oo l 

D[u>]=- I drtr R i SR - — cxp(- i Q L Qt) L ed [uj] Q p TE cxp(z uj t) 

Jo 1 — 2j 

r°° r i 

-i dr dstr K L SR - — exp(- i Q L Q r) Q L SK Vexp{i L (r - s)}L ed [uj] p TE exp(z lu s), (88) 
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where S is given by 

poo 

S = S(oo) = 1 -0 _1 (oo) = -i / rfrexp(-iQLQr) QL SR Pexp(iLr). (89) 

Jo 

The summation in (87) is over the same frequencies as those in (8). By substituting the formal solution of Eq. 
(85): 

pi[w] - (i(L a -u)-C)~ 1 {-iL ai [u]p + D[u]} (90) 

into (67) and by using the definition (66) of the admittance for the interaction (8) of the system with the external 
driving field, which is turned on adiabatically at initial time t = 0, we can obtain 

X- LE M - tr M i{ ^ u) _ c {{ \A, , PoHDfr]}, (91) 

which is a general form of the admittance obtained using the TCLE method, where the interference term Dj[w] is 
given by 

1 f°° 1 
Dj [lo] = - / dr tr R L SR 1 ^ exp(- i Q L Q t) Q [ Aj , p TE ] exp(i lo r) 

i f°° f T 1 

+ -/ dr / dstr R i SR - — exp(- i Q L Q r) Q L SR P exp{i L (r - s)} [ A 3 ■, p TE ] exp(i lo s). (92) 

"Jo JO L — X, 

The first term of Dj[lo] represents the effects of initial correlation of the system and reservoir, and its second term 
represents the memory effects. 

Expanding Eq. (85) up to second order in powers of the system-reservoir interaction 7Y SR , it reduces to 

- i lo Pl [w] = - i L s Pl [lo] + C (2) pi [lo] - i L ei [lo] (p s + p 2) ) + [lo] , (93) 

where the collision operator and the interference term D^ 2 ) [lo] are given by 

/>oo />oo 

C'W = C < 2 ) (oo) = - / dr (L SR exp(- i L r) L SR cx P (z i r)) R = - / dr (L SK i SR (- t)) r , (94) 
Jo Jo 

pOO pT 

i> (2) [lo] — L> (2) [lo] - i j dr ds tr R L SR exp(- i i t) L sr cxp{i L (t - s)} L ed [uj] p s p R exp(i lo s), (95) 
Jo Jo 

with Z)( 2 )[w] given by (74) or (75). The first term of D^[lo], which is called the "first interference term", is equal 
to the interference term [lo] in the TC equation and represents the effects of initial correlation of the system and 
reservoir. Its second term, which is called the "second interference term" , represents the memory effects. Substituting 
the formal solution of Eq. (93) : 

pM = ( i (i s - w) - y'i-i iedH (p s + P^) + £ (2) H } (96) 
into (67) and using the definition (66) of the admittance, the admittance takes the form [41, 43] 

CM - — J } _ g(2) {i [Aj , ps + ] + D?M), (97) 

which is a general form of the admittance obtained using the TCLE method in the lowest Born approximation for 
the system-reservoir interaction. Here, the interference term £r. [lo] is given by 

D^[lo] = D^[lo] + - dr I ds tr R L SR cxp(- i L Q r) L SR exp{i i (r - s)}[A, , p s ] p R exp(iws), (98) 
h Jo Jo 

— (2) — (2) 

with £r. given by (47) or (78). The first term of D y - '[lo] is equal to D y - '[lo] and represents the effects of initial 
correlation of the system and reservoir. Its second term represents the memory effects. The above result (97) can be 
obtained by expanding the TCL equation (82) up to second order in powers of the system-reservoir interaction too 

f A t A Ol 
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IV. RELATIONS AMONG ADMITTANCES 

In this section, we examine the relations among the forms of admittance derived in Sections 2 and 3 analytically. 
As found in the previous section, the admittance (38) obtained using the relaxation TC method has the same form 
as the admittance (72) obtained using the TCE method, i.e., 

x?J c M = xIfM. (99) 

We examine analytically the relations among the admittances (44) [or (77)], (60) and (97), which are obtained using 
the relaxation TC method [or the TCE method] , the relaxation TCL method and the TCLE method, respectively, in 
the lowest Born approximation for the system-reservoir interaction. 

One of the authors has investigated the relation between the admittance (60) obtained using the relaxation TCL 
method and the admittance (97) obtained using the TCLE method [41-43]. We first survey the relation between the 
admittances (60) and (97) in a compact way. The first term of the admittance XiJ CL (<^) [(60)] can be rewritten as 

pOO pt 

(first term of X*f L M [(60)]) = j- J dttr A* exp^ji (w - L s ) t + t + J dr (C™ ( T ) - ) } [Aj , p s + p™ ] , 

= ^J q dttrAiU(u,t)exp_{J dr U~\u}, r) (C (2) (r) - C (2) ) U(u, t)} 



[Aj, p s +p 2) ], (100) 



with U(u,t) defined by 



U{u, t) = exp{ i(u-L 3 )t + C {2) t }, (101) 
and can be expanded in powers of the system-reservoir interaction as 

poo 

(first term of x*J CL M [(60)]) = - J dt tr A, U(u, t) [A, , p s + p (2) } 

p OO pt p oo 

/ dt dr dstvAi U(u,t — t) tr R i SR exp(— i L s)L S n.ex.-p(i L s) PhU(uj,t) 

JO Jo Jt 



I 

h 



x [Aj ,p s + p M ] + ..., (102) 



where (58) and (94) have been substituted into (t) and , respectively. The first term of (102) can be integrated 
to give 

poo 

- J dt tr Ai U(uj, t) [Aj , p s + p {2) ] = - tr At ( i (L s - w) - )~ l [Aj , Ps + p {2) ] . (103) 
Considering that in the second term of (102), 

pOO pt pOO poo poo p oo 

/ dt dr dstr AiU(uj,t - t) = / dr dt dstr AiU(u,t - r), 

JO JO Jt Jo Jt Jt 

pOO pOO pOO pS 

= dr dstrA^iiLs-^-C^y 1 = / ds dr tr A, ( % (L s - u) - C (2) ) ~\ (104) 
Jo Jt Jo Jo 

the first term of the admittance x^J CL (o>) [(60)]) can be expressed as 

(first term of x*J CL M [(60)]) =trA t (i(L s -u>)- C {2) ) _1 {{i/K)[Aj , p s + p (2) } + (second term of D (2) [lo])} 

+ HOT (higher order terms), (105) 

(2) 

where D) '[uj] is given by (98). The second term of the admittance (60) can be rewritten by substituting (59), as 
(second term of x*f L M [Q0}) = -^J dt J rfrtr A; exp^ { (i (w - L s ) + C™)(t - t) + J ds (C (2) (a) - C (2) )} 

x / d(3' tr K L S Kexp(-iL T)[Aj , p s fh.U s ^{- ihf3')] exp(iwr), (106) 
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and can be expanded in powers of the system-reservoir interaction as 

-i />oo poo pf3 

(second term of x* TCL H [60]) = - ^ J dr J dt J d/3' tr A, cxp {(* (w - L s ) + C (2) ) (t - r)} 

x tr Rj L SR cxp(-iL r) [Aj , p s p R H SR (- i h(3')] exp(iwr) 
+ ■■■, (107) 

which can be expressed by integrating with respect to t, as 

(second term of x* TCL M [(60)]) = tr A l (i(L s -w)- C* (2) ) _1 (Df } [w] or first term of Df ] {uj]) + HOT, (108) 

where Z)j 2 ^[w] is given by (47) and coincides with the first term of D^ 2 \u\. Expansions (105) and (108) show that 
the admittance obtained using the relaxation TCL method, coincides with the admittance obtained using the TCLE 
method in the lowest Born approximation for the system-reservoir interaction [42, 43], i.e., 

x*f L H [(60)] = xI, CLE M t( 97 )] + H0T ( hi s hcr order terms )- ( 109 ) 

We next examine the relation between the admittance (60) obtained using the relaxation TCL method and the 
admittance (44) obtained using the relaxation TC method [or the admittance (77) obtained using the TCE method]. 
The collision operator [u>] given by (46) can be rewritten as 

poo 

C {2) [u] = - dTtr R L SR cxp(-iLoT)L SK p R cxp(iLJT) = C {2) +AC (2) M, (110) 
Jo 

with given by (94), where AC* 2 '[w] is defined by 

/•OO 

AC (2) [u] = - / dTtT R L SR cxp(-iL T)L SR {exp(iujT) - exp(z L t)} p R . (Ill) 
Jo 

Then, the first term of the admittance XiJ CL (w) given by (60) can be rewritten as 
(first term of X *f L M [(60)]), 

POO /■*_-! 

= -J dttrA l exp^[t(cu-L s )t + C^[u}t + J dr(C( 2 )(T)-C^H)}[ J 4,, /9s + )0 [ ) 2) ], 
= l -J dtt r A l U(u J ,t)exp^{J drU-\u J ,T)(C (2 \T)-C^-AC^\iu})U(u J ,T)j[A J ,p s + p { 2) ], (112) 
with U(u),t) defined by 

U((j, t) = exp {i {u - L s ) t + C {2) [u] f } , (113) 
and can be expanded in powers of the system-reservoir interaction as 

p OO 

(first term of xff'M [(60)]) = - J dt tr A l U(w, t) [A, , p s + p (2) ] 

dr tr A, U{lj, t - r) { J ds tr R L SR exp(- i L s) L SR cxp(i L s) p R 
+ J dstr R i SR exp(- iL s) L SR {cxp(i ws) - exp(zL s)} p R | 



+ 



xU{uo 1 T)[A J , Ps +p (2) } + ■■■ , (114) 



where (58), (94) and (46) have been substituted into C (2) (t), and C (2) [w], respectively. The first term of (114) 
can be integrated to give 

/OO 
dt tr Ai U(w, t) \Ai , oa + i 2) ] = - tr A: (i(L«-u>)- (7 (2) U V 1 \A< , p s + 1 . (115) 
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The second term of (114) becomes, by integrating in the same procedure as in (104), as follows 

f°° l ( f r 

J dr tr Ai ( i (L s - ui) - (7 (2) [u] ) |y ds tr R L SR cxp(- % L t) L sr cxp(z L r) U(u), s) p R 

+ J ds tr R L SR cxp(- i L r) L SR {cxp(iwr) - exp(iL r)} E/(w,s) p R j [Aj , p s + p 2) ] 

■ / dr tr ^ ( i (L s - w) - C {2) [u] ) tr R L SR exp(- i L r) L SR { exp(i L t) exp ( i (u> - L s ) r + C (2) [w] r ) 
Jo 

- cxp(i L r) - (cxp(z w r) - cxp(z L r)) } ( i (w - L s ) + C* (2) [w] ) p R [A,- , p s + Pq 2) ] > ( 116 ) 



i 

ft 



i 
ft 



which is the higher order term in powers of the system-reservoir interaction. Thus, the first term of the admittance 
X^J CL (w) [(60)] can be expressed as 

(first term of x* TCL H [(60)]) = tr A, (i (L s - w) - C( s '[w]) _1 (i/ft) [A,- , p s + p[, 2) ] + HOT. (117) 
The second term of the admittance XiJ CL (w) given by (60) can be rewritten by substituting (59), as 
(second term of x* TCL H [(60)]), 

= ~K jf ^ /' dT jf d/3 ' tr ^ CXP ^ { ( z ( w " L s) + C {2) [u] ) (t - t) + J* ds (C (2) (s) - <7 (2) [w]) } 

x tr R L SR exp(- i L r) [Aj , p s PrWsr(- i ft/3') ] exp(zwr), (118) 
and can be expanded in powers of the system-reservoir interaction as 

i r-oo roc r/3 

(second term of Xy^M [60]) = -- / dr dt d(3' tr A l exp {(i (w - L s ) + C (2) [w] ) (t - r) } 

" JO Jr JO 

x tr R L SR exp(-iL T)[Aj , p s Pr H SR (- « ft/3') ] exp(iwr) H , (119) 

which can be expressed by integrating with respect to t, as 

(second term of x*/ CL M [(60)]) = tr A (i (L s - u) - C^w]) -1 Df^cj] + HOT, (120) 

with D^p[uj] given by (47). Expansions (117) and (120) show that the admittance obtained using the relaxation TCL 
method, coincides with the admittance obtained using the relaxation TC method [or the admittance obtained using 
the TCE method] in the lowest Born approximation for the system- reservoir interaction, i.e., 

X*J CL M [(60)] = x?f (w) [(44)] (or X Jf M [(77)]) + HOT (higher order terms). (121) 

We also examine the relation between the admittance (44) obtained using the relaxation TC method [or the 
admittance (77) obtained using the TCE method] and the admittance (97) obtained using the TCLE method. The 
operator (i (L s — uj) — C^[u]) can be rewritten by using AC' 2 '[w] given by (111), as 

(i (L s - lo) - C^uj])- 1 = {i (L s - u) - - AC< 2 >M}-\ 
= {(i (L s - oj) - C^) {1 + (i (L s - w) - C^)- 1 (- AC (2) H)}} _1 , 

= {1 + (i (L s - w) - C^)- 1 (- AC^M)}" 1 (i (L s - w) - C (2) )~\ (122) 
and can be expanded in powers of the system-reservoir interaction as 

^(Ls-^-C^l^y 1 = (t(L s -u;)-C^y 1 {l + AC^[u J }(i(L s -Lu)-C^y 1 }+ ••• . (123) 
Substituting the above expansion into x^J c (w) given by (44) [or xJj E {u) given by (77)], we have 

X »J C M =trA l (i (L s - u) - C^y 1 {(i/h) [Aj , p s + } + Df[u]} 

-+- tr 4, d(L, - - n^A' 1 An ( V\,,,] d(L, - - aW)' 1 u.ik\ \a, . n, + «1 2) 1 + hot. c\ia\ 
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where is given by (47) or (78), and is equal to the first term of the interference term Z?j 2 ^[w] given by (98). 

Substituting (111) into the second term of the above expansion and considering that 

AC^[w] (i (L s -w)- C^Y 1 {i/h) [Aj , p s + /4 2) ] , 
= - rfrtr Rj L S RCxp(-zL r) L SK {exp(swr) - exp(iL i~)} p R (i (L s - uj) - C (2) ) 1 - [Aj , p s + p{, 2) ] , 

= T dr ds tr a L SR exp(- i L r) L SR exp{« L (r - s)} p R [Aj , p s ] exp(i u> s) + HOT, 
n- Jo Jo 

= (second term of Df ] [u]) + HOT, (125) 

we find that the admittance obtained using the relaxation TC method [or the admittance obtained using the TCE 
method], coincides with the admittance obtained using the TCLE method in the lowest Born approximation for the 
system- reservoir interaction, i.e., 

X*f M [(44)] = X £V) [(77)] = x!f E (u,) [(97)] + HOT (higher order terms). (126) 

According to the relations (99), (109), (121) and (126), the admittance obtained using the relaxation TC method 
[or the TCE method] , the admittance obtained using the relaxation TCL method and the admittance obtained using 
the TCLE method, coincide with each other in the lowest Born approximation for the system-reservoir interaction 
Hsu, i.e., 

x£V) = xgV) - x*f L M - XSTM, ( U P to °«)) ( 127 ) 

though the exact admittances (38) [or (72)], (56) and (91), which are obtained using the relaxation TC method [or the 
TCE method], the relaxation TCL method and the TCLE method, respectively, take the forms different from each 
other. The admittances that are obtained using the relaxation TC method [or the TCE method] , the relaxation TCL 
method and the TCLE method, respectively, in the lowest Born approximation for the system-reservoir interaction, 
have the same second order term and the higher order terms different from each other. 

V. HIGHER ORDER EXPANSIONS 

In this section, we derive the formulae necessary for the higher order expansions in powers of the system-reservoir 
interaction Hsr, and give the forms of the admittances in the n-th order approximation for Hsr- 

A. Expansions of the collision operators C(t), C and C[u>] 
The collision operator C(t) given by (53) can be rewritten as 

oc 

C(t) = i (L SR £(t) {1 - X(t)}-\ = ~ * E ( Ls * E (*)">b > ( 128 ) 

71=1 

with T,(t) defined by (84). Defining the unperturbed evolution operator Uo(t) as 

Uo(t) = exp(- i L 1) = cxp{- i (i s + L K ) t}, (129) 

we have 

L SR (-t) = U (t)L SR U-\t), H SR (t) = (UoWHm) =cxp(» Hot/h) W SR cxp(- z H t/h), (130) 

and 

exp(iLt) = I/^ 1 (t)exp < _{* J drL SR (-r)}, (131) 
exp(-iQiQt)Q = exp^[-i f dr Q L SR (- r) q} U (t) Q . (132) 
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Then, £(£) defined by (84) can be rewritten as 

£(i) = -ijf drexp^j-i^ ds 1 QL SR (-.s 1 )Q}QL SR (-r)7'cxp^{i^ ds 2 L SR (- s 2 )}, (133) 

which can be expanded in powers of the system-reservoir interaction Hsr as 

E(t) = (i) + £< 2 ) (t) + £( 3 > (£) + •••, (134) 

where Y,( n \t) is the n-th order part in powers of Hsr- The first-order, second-order and third-order parts of E(£) take 
the forms, respectively, 

W(t) = -i f drQL SK (-T)V, (135) 
Jo 

^{t) = {- if f dn f 1 dr 2 {QL SE (- r 2 ) Q L SR (- n) P - Q i SR (- n) ? L SR (- r 2 )}, (136) 
Jo Jo 

Z (3 \t) = (-i) 3 f dn [ 1 dr 2 f 2 dr 3 { QL sr (-t 3 )Q J L SR (-r 2 )QL SR (-r 1 )7' + QL SR (-T 1 )P J L SR (-r 2 )L SR (-r3) 
Jo Jo Jo 

- Q Lse(- r 2 ) Q L SR (- n) P L SR (- r 3 ) - Q L SR (- r 3 ) Q L sa (- n) 7? L SR (- r 2 )} . (137) 
In the derivation of T,^(t), the following integral transformation has been used 

/■Ti /-Ti /-Ti /-T2 /"Ti /-Ti /'Tl 1>T 2 rTl 1>T 3 

/ c?r 2 / dr 3 = dr 2 / dr 3 + / dr 2 / dr 3 = dr 2 / dr 3 + / c?t 3 / dr 2 . (138) 

JO JO JO JO JO Jr 2 JO JO JO JO 

Expanding the operator C(t) in powers of the system-reservoir interaction Hsr as 

C(t) = C (2) (t) + C {3) (t) + C {4) (*) + •••, (139) 

where (t) is the n-th order part in powers of Hsr, the second-order part of C(t) is given by (58) and the third-order 
and fourth-order parts can be obtained as 

(t) = (- if J dn J* dr 2 (L S r Lsr(- t 2 ) Lsh(- n)) R , (140) 

C {i \t) = (-z) 4 f dn [ Tl dT 2 P dr 3 {(LsRLsR(-T 3 )LsR(-T 2 )LsR(-n)) R -(LsRLsR(-T 3 )) R (L S R(-T 2 )LsR(-n)) R 
Jo Jo Jo 

- (LsrL S r(-t 2 )) r (L S r(- T 3 )L SR (-n)) R - (L S r L S r{- ti)) r (L S r(- t 3 )L sr (- t 2 )) r }, 

(141) 

where the following integral transformation has been performed 

ft pT\ ft ft fT\ [-T2 PT\ ft 



dn I dr 2 / dr = dn / dr 2 < dr + dr + dr>, 
Jo Jo Jo Jo *■ Jo Jt 2 Jti - 1 

/■t /-Tl /-T2 ft f-T\ r-T ft r-T r-T\ 

= dn dr 2 / dr+ dn dr dr 2 + dr dn dr 2 . (142) 
Jo Jo Jo Jo Jo Jo Jo JO Jo 



In the derivations of (140) and (141), we have considered that (Hsr(— i)) R = by virtue of the rcnormalization (7). 
The n-th order part C^ n \t) in the expansion of the collision operator C(t) has the following structure: 

C< n >(t) = (-*)" / dn I " dr 2 f 2 dr 3 --- f " 2 dT„_ 1 ( J L SR (0)L SR (-r„_ 1 )L SR (-r„_ 2 )...i SR (-T 2 )i SR (-T 1 )) oc , 
Jo Jo Jo Jo 

(143) 

for n > 2, where (• • • ) oc denotes the "ordered cumulants" [48-50]. 

The expansion formulae of the collision operator C can be obtained by expanding C [= C(oo)] in powers of Hsr as 
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The collision operator C[u>] given by (34) can be rewritten by considering that (W SR ) R = 0, as 

poo ft 

C[w] = - dt(L SK cxp^{-i dTQL SR (-T)Q}U (t)L SR ) R exp(iwi). (145) 
Jo Jo 

Expanding the collision operator C[lo] in powers of the system- reservoir interaction Hsr as 

C[u] = C (2) H + [u] + [u] + ■ ■ ■ , (146) 

where [cj] is the n-th order part in powers of Hsr, the second-order part of C[u>] is given by (46) and the third-order 
and fourth-order parts can be obtained as 

/•CO nTl 

CW[u] = (-i) 3 J dnj dr 2 (isRisR(-r 2 )L SR (-T 1 )) R cxp{ i (^-L s )r 1 }, (147) 

r CO P Tl P T2 

C (4) H - (- i) 4 / dn / dr 2 / dr 3 {(L S r L SR (- r 3 ) Lsa(- r 2 ) L SR (- n)) R 
Jo Jo Jo 

- (L S r Lsa(- t 3 )) r (L sr (- r 2 ) L SR (- n)) R } exp{i (w - L s ) n} ■ (148) 
The n-th order part C'^^w] in the expansion of the collision operator C[uj] has the following structure : 

C^M = (-»)" /*dTi r dr 2 f T ' 'dr 3 --- f" 2 dr„_ 1 <L SR (0)L SR (-r„_ 1 )i SR (-r„_ 2 )...L SR (-r 2 )i SR (-r 1 )) 
Jo Jo Jo Jo 

x exp{i (w - L s )ti} , (149) 
for n > 2, where (• • • ) pc denotes the "partial cumulants" [48-50]. 

B. Expansions of the inhomogeneous terms and interference terms 

We first expand the thermal equilibrium state Pte given by (3) in powers of the system-reservoir interaction Hsr- 
Writing the expansion of exp(— (3H) = cxp{— (3 (Ho + Hsr)} in powers of Hsr as 



exp(- P H) = exp(- Ho) exp^ { - / d(3 ' exp(J3 ' Ho) H SK exp(- j3 ' Ho) } = exp(- p Ho) { 1 + ]T S„ } , 

"' n=l 



(150) 



with ^„ dehned by 

S„ = (-)" /^V f 1 d8 2 --- f n 1 dpnHszi-ih^Hszi-ihfo) ■■■H SR (-ihp n ), (151) 
Jo Jo Jo 

we have 

Trexp(-/3W) = {Trexp(-/3Wo)}-(l + <e) ; £ = &+&+& + •■■, (152) 

where we have defined 

£„ = Tr{cxp(-/3H )S Il }/(Trcxp(-/3H )) - Trp s/0R £„ , 

= (-)" / dfr f 1 dft • • • / " 1 d/3„ < Wsb(- i fift) Wsb(- * Tift) • • • Wsb(- j , (153) 

Jo Jo Jo 

with the notation (• • • ) SR defined by (• • • ) SR = Tr p s p R ■ • • , where p s and p R are given by (4) and (5), respectively. Here, 
we have considered that (Hsr(— ifi/3)) R = 0. Then, p TE can be expanded as 

Pte = exp(- /3 Ho) (1 + Hi + S 2 + S 3 + • • • ) / {( Tr cxp(- /3 W„))-(l + 0}, 

= p s p R (1 + E 1 + E 2 + E 3 + ■ ■ ■ ) (1 - £ + - + ■ ■ ■ )• (154) 

Expanding p TE in powers of the system-reservoir interaction Hsr as 

_(0) , fl) , (2) , (3) , /1rr x 
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where is the n-th order part in powers of Wsr, we have 



rP 

Pte = Ps Pr , Pte = Ps Pr Si = - p s Pr / d/3i H SR {- ih/3i), (156) 

Jo 

rP rPi 

Pte=PsPr(S2-&) = PsPr / dft / dp 2 {H SR (-ihp 1 )H SR (-ihp 2 ) - (7M-*W«sr(-*W) sr }, (157) 

Jo Jo 

Pte = Ps Pr (S 3 - 6 ~ Si £ 2 ), 

= -PsPr f d0! f'dfo [' df3 3 {n SK (-thf3 1 )H SK (-thf3 2 )H SK (-th(3 3 ) 
Jo Jo Jo 

- (W S r(- * ft/3i) W S r(- t ftft) Wsr(- i } 

+ pspr / dp'n SR (-ihp') d/h dp 2 (n SR {-ihp 1 )n SR {-thp 2 )) SK , (158) 

Jo Jo Jo 

Pre = Ps Pr { S 4 - 6t - (S 2 - 6) 6 - Si £ 3 }, 

= PsPr /"d/3i /^dft /^dft dp 4 {n SR (-ihp 1 )n SR (-ihp 2 )n SR (-ihp 3 )n SR (-ihp 4 ) 
Jo Jo Jo Jo 

- ihfr) H SR (- * ft ft) Wsr(- » ftft) Wsr(- * W) sr } 

/■/? /-ft 

- Ps Pr / dfc / d/3 2 {W SR (- * hpi) H SR {- i /ift) - (W S r(- ihpi) H S r(- ihp 2 )) sz } 

Jo Jo 

rP rP' 

x dp' dp"(H SR (-ihp')H SR {-ihp"))^ 
Jo Jo 

rP rP rPi rPi 

- PsPr J dp"H SR (-ihp') J dp, J dp 2 J^ dp 3 {H SR {-ihp 1 )Hs R {-ihp 2 )Hs R {-ihp 3 ))^. (159) 

The inhomogeneous terms I(t) and I(t) which are given by (31) and (54), can be rewritten using (132) as 

I(t) - -i-tr R L SR exp^{-i drQi SR (-r)Q};7o(t)Q[A, pte], (160) 

I(t) = - i-tr B L SR {1 - S(t)}- 1 exp^j - i jf dr Q L SR (- r) q} C/ (t) Q [A , p TE ] , (161) 

which can be expanded in powers of the system-reservoir interaction H SR , respectively, as 

J(t) = 7 (2) (t) + 7 (3) (t) + 7 (4) (*) + •■■, J(t) = / (2) (t) + / (3) (t) + / (4) (*) + ■■■, (162) 

because /W (i) = 1^ (t) = according to Q p$) = Q Ps Pr = 0, where and /(")(*) are the n-th order parts in 

powers of Hsr- The second-order parts are given by (42) and (59), and the third-order and fourth-order parts can be 
obtained as 

,t ,/3 



p)(t)=j(3)(t)= / dr / d/3itr R L SR L SR (-T)[/oM[^ PsPrW S r(-^/?i)] 
Jo Jo 

/■/3 /■/?! 

-2/ d/3i / d/3 2 trRL SR C/ (i)[A, Ps/9 RH SR (- t ft/3i)H S R(-ift/3 2 )], (163) 
Jo Jo 

/ (4) W=-* / dn [ 1 dr 2 [ P dp 1 tr R L S RL S R(-T 2 )QL S R(-T 1 )U (t)[A, pspRHsRi-ihP!)] 
Jo Jo Jo 

ft rP rPi 

- dr dP\ I dp 2 tr R L SR L SR (-T) 
Jo Jo Jo 



x U (t) [A , Ps p R {H SR (- i hp,) H S r(- ihp 2 )- (H S r(- ihfa) H SK (- ihp 2 )) K }} 



+ 1 f dp, f 1 dp 2 f 2 dp 3 tT K L SK U a (t)[A, psPRH^-ihp^H^-ihp^H^i-ihps)] 
Jo Jo Jo 
rP rP fPi 

-i dp' dp, dp 2 tr R L S RU (t)[A, p sPR H S R(-ihp')](H S R(-ihp 1 )H S R(- l hp 2 )) SR , (164) 
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lW(t) = n 4 \t)+i [ dn [ 1 dr 2 f ' dp 1 tr R L SR L SR (-T 1 )TL SR (-T 2 )U (t)[A, p fl p tL 'H m (-ihfa)], (165) 
Jo Jo Jo 

where we have considered that (Hsr(— t)) R = (Hs R (— iftP)) R = (Hsk)r = and Q p^) = p[ E . The second term of (165) 
comes from E(i) = 1 — 6 (t) which produces by renormalizing the time-convolution and represents the memory effect. 
Thus, the fourth-order part I^ 4 \t) of I(t) includes the memory effect. 

The interference terms D[lo] and D[uj] which are given by (71) and (88), can be rewritten using (131) and (132) as 
D[uj}=-J^ drtr R L SR cxp^{-z^ ds Q L SR (- s) q} U (t) L ed [uj] Q p TE e iu T , (166) 

p OC p T 

D[w}=- drtr R L SR (^E")exp^{-i / ds QL SR (- s) q} U q (t) L ed [uj] Qp TE e luJT , 
Jo n=0 Jo 

-ijf dnjf dr 2 tr R I, SR (^E")exp^{-^ dsi Q L SR (- si) q} U (n) Q L sr 



cxp^ji / ds 2 L SK (- S 2 )j LedM PTEGVp(iuJT 2 ), (167) 

which can be expanded in powers of the system-reservoir interaction 7Y SR , respectively, as 

D[u] = £ (2) [ u ] + £»(3) [ w ] + flW [ w ] + . . . , D[u] = [lu] + [lu] + D< 4 ' [u] + ■ ■ ■ , (168) 

because [to] = [lo] — according to Q p$) = Q p s p R = 0, where [lo] and [lo] arc the n-th order parts in 
powers of Hsr- Here, E is given by (89) and can be expanded in powers of H SR as E = E(oo) = E^ + E^ 2 ) + E< 3 > + • • • 
with SW = E(™)(oo). The second-order parts are given by (74) [or (75)] and (95), and the third-order and fourth-order 
parts can be obtained as 

/•OO PT\ 



/•OO /■Tl pp 

D ( - 3) [lo]=-i / dn / dr 2 / dpi tr R L SR L SH (- r 2 ) f7o(n) £ e dM Ps Pr H S r(- i ftPi) exp(i w n) 
Jo Jo Jo 

/•oo /./3 /.ft 

/ dr / d/?i / d/3 2 tr R L SR C/ (t) £ e d[^] p s Pr Wsr(- « ft/?i) W SR (- i ft/? 2 ) exp(i wr), (169) 
Jo Jo Jo 

/•OC /■Tl 

D (3) [w] =D (3) [u] - / dn / dr 2 / dr 3 tr R L SR {L S r(- t 3 ) L sr (- n) Uq(t 2 ) L b6 [lo] p s p R exp(i lo t 2 ) 
Jo Jo Jo 

+ L sr (-t 2 ) L SR (- ti) J7 (r 3 ) £ ed [w] ps Pr cxp(i w r 3 )} , (170) 

/'CO /-Tl /-T2 

D (4:) [lo]=- dn dr 2 dr 3 / dpi tr R L SR £ S r(- t 3 ) Q L SR {- r 2 ) C/o(n) £ ed [w] p s Pr H SR (-i ft/?i) exp(i wn) 
Jo Jo Jo Jo 

+ i / dn / dr 2 / d/3i / d/3 2 tr R L SR L SH (- r 2 ) C/ (n)i 8 dM p s PRexp(i w n) 
Jo Jo Jo Jo 

x {H SR (- z ft A) H m (- 1 ft/3 2 ) - (Wsr(- * ft/3i) W SR (- i ft/3 2 )) R } 

/■oo /■/3 /■/3i /■/3 2 

+ / dr / d/3i / d/3 2 / d/3 3 tr R L SR [/ (r)L ed [w] p s Pr H S r(-« ft/3i)W SR (-i ftP2)U S R.{-i Hfo) exp(i w r) 
Jo Jo Jo Jo 

/■oo /-/3 /-/3 /■/3i 

-/ dr dp' dPi dp 2 tr R L SR Uo(T)L ed [u J }p s p R H S R(-iftp r )(H S R{~iftp 1 )H SR {-iftp 2 )) 
Jo Jo Jo Jo 



SR 

x exnf?' i,} rV fl 71 1 
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D (4) [w] =5 (4) [w] 

/>00 (-00 rT\ r/3 

+ dr dn dr 2 / dft tr n L sn L SR (- n) 7> £ SH (- r 2 ) f/ (r) £ e dM p s p R W SR (- i ft ft) exp(i wt) 
Jo Jo Jo Jo 

/>0O /■T 1 /"T 1 Z* 5 ! 

+ i / dn / dr 2 / dsi / ds 2 tr R L SR L SR (- s 2 ) QL SR (- L SR (- n)U a (T 2 ) L ed [uj} p s p R cxp(iujT 2 ) 
Jo Jo Jo Jo 

/>0O /-Ti />OG />Si 

-i/ dn dr 2 dsi I ds 2 tr R L SR L SK (- si) V L SR (- s 2 ) L SR (- n) U (t 2 ) L ed [ui] p s p R exp(i wr 2 ) 
Jo Jo Jo Jo 

/•OO /-Tl rTl—T-2 PT3 

+ i dn dr 2 dr 3 dntr R L SR U Q (n) L SR V U^ 1 ^ - t 2 ) L SR (- t 3 ) L SR (- n) L ed [u;} p s p R 

Jo Jo Jo Jo 

x exp(zwT 2 ) 

/■oo /■ri PT1—T2 rf) 

l dn dr 2 dr 3 dft tr R L SR f/ (n) ^sr 7> {/^(n ~ r 2) ^sr(- 73) L ed [uj] p s p R H SR (- i ftft) 

Jo Jo Jo Jo 

x cxp(i cj t 2 ) 

/•oo pri p(3 pfSi 

-i dn dr 2 / dft / dft tr R L SR L SR (— n)VUo(T 2 ) L ed [u)\ Ps Pr cxp(i uj t 2 ) 
Jo Jo Jo Jo 

x {H S R(-iftft)W SR (-iftft) - (H SR (-ihf3 1 )H SR (-ih(3 2 )) SR }, (172) 

where we have considered that (Hsr(— t)) R = (Hs R (—ifrP)) R = (Wsr)r = and Q p[% — p^ . The interference term D[lu] 
consists of the interference term D[u>] in the TC equation, which represents the effects of initial correlation of the 
system and reservoir, and of the terms that represent the memory effects. 

C. Admittance in the ri-th order approximation for 7^sr 

In this subsection, we give the forms of the admittances in the n-th order approximation for the system-reservoir 
interaction Hsr- The admittance obtained using the relaxation TCL method by expanding Eq. (52) up to n-th order 
in powers of Hs R and by proceeding in the same way as in Subsection III.B, takes the form 

XT L (")=1 dttvA i exp^{i(u;-L s )t+ dr YJ C^\r)} A,- , Ps + YJ p^ ] 

^° ^° m=2 m=2 

■ roc /•* -t n n 

+ -/ dt dTtrA i exp < _{i(w-L s )(i-r)+ / ds C (m) (s)\ YJ /j m) (r) exp(iw t), (173) 
h Jo Jo [ Jr m=2 J m=2 



with p m ' = tr R p^ , where /j™' (t) is equal to (t) with Aj in place of A. The admittance obtained using the 
relaxation TC method by expanding Eq. (28) or (32) up to n-th order in powers of Ti SR and by proceeding in the 
same way as in Subsection III. A, is equal to that obtained using the TCE method by expanding Eq. (63) or (68) up 
to n-th order in powers of Hs R and by proceeding in the same way as in Subsection III.C, and takes the form 

n _^ . n n 

x£V) = xlf M = trA i ^ - «) - E ° {m) ^) " {\ • (>* + E p ( o m) } + E £} m) M} ■ ( 174 ) 

m— 2 m— 2 m— 2 

The second-order part Z)j 2 ^[u;] of is given by (47) or (78), and the third-order and fourth-order parts are given 

by 

-Oj 3) M=T / dn / dr 2 / dpitT R L SK L SR (-T 2 )U (Ti)[Aj , p s p R HsK{-ihPi)]exp(iu>n) 
11 Jo Jo Jo 

+ -/ dr / dft / dft tr R i S R^o(r)[^-,PsPRWsR(-i ftft) W SR (-i ftft)] exp(iwr), (175) 
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l Z" 00 /" Tl Z" 1 " 2 f 13 

Dy[u)]=-j dn J dr 2 J dr 3 J d/3i{tr R L SR i SR (- t 3 ) L SR (- r 2 ) - (i SR i SR (- t 3 )) r tr R L SR (- r 2 )} 

x U (Ti)[Aj , psPRW SR (-ifi/?i)] exp(iwTi) 

~ T. I dTl dT2 / d ^ x I dp 2 tT R L SR L SR (-T 2 )U (Ti)exp(iuTi) 
" Jo Jo Jo Jo 

x , p s p R {H SR (- ihfr) Hsn(- ihfo) - (Wsa(- iSft) W SR (- ^/? 2 )) R }] 

1 /-oo /-/3 r&\ rfh 

--/ dr dfr df3 2 d(3 3 tT R L SR U (T)[A ]7 pspiHto{-ihp 1 )H s &(-ihp 2 )H WL {-ihfc)] 
n Jo Jo Jo Jo 

x exp(iwT) 

+ hJo ^ Jo Jo dPl Jo d fc tr * Ls * U °( T ^ A ^ P^^sK(-ihf3 r )](H SR (-ih(3 1 )n SR (~zhf3 2 )) SR 

xcxpfiwr). (176) 



The admittance obtained using the TCLE method by expanding Eq. (82) or (85) up to n-th order in powers of Hs R 
and by proceeding in the same way as in Subsection III.D, takes the form 

n _^ . n n 

X ™(u) = tvA, (i (L S - U ) - £ C {m) Y {~ h h Po n) ] + £ ^ m) H} ■ (177) 

m— 2 m— 2 m— 2 

The second-order part D^[u;] of Dj[uj] is given by (98), and the third-order and fourth-order parts are given by 

£>f ) M = Df\u] + - / dri / dr 2 / dr 3 {<L sr L sr (-t 3 )L sr (-t 1 )) r [A,- (-7*), Ps ] exp(iwr 2 ) 
ft Jo Jo Jo 

sr ^sr(— t 2 ) L SR (— Ti)) R [A j (-T 3 ),p s ]exp(iwr 3 )}, (178) 



r-oo /.oo 



ft. 



/-00 /-00 /.Tl /.p 

- / dr / dn / dr 2 / (i S RisR(-Ti)) R tr R i SR (-r 2 ) U (t) [A 3 ■, p s p R H SR (- i fr 0i)] exp(iwr) 
n Jo Jo Jo Jo 

£ I-OO /-Tl /-Tl /-T3 

- (L SR L SR (- t 4 )) r (L sr (- t 3 ) L SK (- ^^[Aji- t 2 ) , ps] exp(iujT 2 ) 

p OO pT\ p OO p T3 

■ / dri / dr 2 / d-r 3 / dr 4 (£ S r £sb(- t 3 )) r (L S r(- t 4 ) £ S r(- n)) R [A,-(- r 2 ) , ps] exp(i w r 2 ) 
Jo Jo Jo Jo 

/•OO /-Tl rTl—T2 rT3 

/ dri / dr 2 / dr 3 / dr 4 (L SR L SR (- n)) R C/ (t 2 ) (£ S r(- t 3 ) L S r(- t 4 )) r [A,- , p s ] exp(i wt 2 ) 
Jo Jo Jo Jo 

+ / dTl / dT2 / dT3 / ^ ( L SR-£'SR(-Ti)) R J7o(T 2 )tr R i SR (-T 3 )[^ J :, /9 S PrW S r(-^/3i)] cxp(jwr 2 ) 
ft Jo Jo Jo Jo 

• p 00 /.ti /./3 /./3i 

+ dn / dr 2 J dfc J df3 2 (Ls R L SR (-T 1 )) R Uo(T 2 ) [A jt ps{<WsR(-ift/3i)WsR(-ifi/92)> R 

- (W m (- * A) Wsr(- * fi/3 2 )) gR } ] cxp(z u r 2 ), (179) 

where A(t) is defined by 

= cxp(« L s t) A = exp{iH s t/h)Aexp(-iH s t/h). (180) 



As expressed using time-integrals alone for the second-order interference terms Z?j 2 ^[a;] and in (78) and 

(98), the inhomogeneous term Ij(t) and the interference terms Dj[u>] and Dj[u>] can be expressed using time- integrals 



22 



(n = 2, 3, 4) of Ij (t) take the forms using time- integrals alone : 

/oo 
dr tr R i SR [A s (- i) , i SR (- r) p s p R ] , (181) 

/j 3) (t) = i dn / rfr 2 tr R L SR L SR (- r 2 ) [A,- (- t) , L SR (- n) ps pn ] 
Jt Jo 

/OO <.T! 
y dT 2 tv K L sn [Aj(-t), L sn (-T 2 )L sn (-T 1 )p s pn], (182) 

I f\ t )=i dT i dT 2 dT 3 {tT K L SK L SR (-T 3 )L SR (-T 2 ) [Aj{-t), Lsr(-ti) PsPr] 
Jt Jo Jo 

- (L SR L SK (-T 3 )) R tr R L SK (-T 2 ) [Aj(-t), L SR (- n) ps Pk]} 

r OO f-T\ r-t 

+ dn dr 2 dr 3 { tr R L SR L SR (- r 3 ) [Aj(-t), L SR (- t 2 ) L SR (- n) p s p R ] 
Jt Jt Jo 

- (L SR L SR (-T 3 )) R tr R [Aj(-t), L SR (- t 2 ) L SR (- n) p s p R ]} 
dn J dr 2 j dr 3 tr R L SR [Aj(-t) , L SR (- t 3 ) L SR (- t 2 ) L SR (- n) ps Pr] 

dn J dr 2 J dT 3 (L SR L SR (-T 2 )) R tT R L SR {-T 3 ) [Aj(-t), L SR (-n)p s p R ] . 



(183) 



The second-order interference terms Dj '[oj] and Dj [u>] given by (78) and (98) can be, respectively, expressed in the 
compact forms 

^ 2) M = _ fe/ dr i / dr 2tr R L SR [Aj(-T 2 ) , L SR (- n) p s p R ] exp(«wr 2 ), (184) 
Jo Jo 

/■oo /-ri 

Df ) [uj] = Df ) [uj] + -J^ dnj dT 2 (L SB L SB (-T 1 )) R [A j (-T 2 ),p B ]exp(iwT 2 ). (185) 

By virtue of the results of Appendixes A and C, the third-order interference term D^- '[ui] given by (175) takes the 
form using time-integrals alone : 

/>0O />T i pT 2 

5j 3) [w] = - - J dn J dr 2 J dr 3 { tr R L SR L SR (- t 3 ) [Aj(- t 2 ) , L SR (- n) Ps Pr ] exp(i u t 2 ) 

+ tr R L SR [A(-t 3 ) , L SR (- t 2 ) L SR (- n) p s p R ] exp(awT 3 ) }. (186) 

The third-order interference term is given by (178) with the above By virtue of the results of 

Appendixes A, C and D, the fourth-order interference terms and D^\oj] given by (171) and (172) take the 

forms using time-integrals alone : 

L>j 4) [w]=^ J dn j dr 2 J dr 3 J dr A { tr R L SR L SR (- r 4 ) £ S r(- r 3 ) [Aj(- t 2 ) , L S r(- n) p s p R ] 

- (isRisR(-T 4 )) R tr Rj L SR (-r 3 ) [Aj(-T 2 ), L SK (- n) Ps Pr] } exp(z u t 2 ) 

i Z" 00 f Tl f~ 2 f Ta 
+ hj dTl J dT2 J dT3 J dT i{ tT ^ L ^ L sR(-T 4 )[Aj(-T 3 ), Ls R (-T 2 )L SR (-n)psPn] 

- (£sr£sr(- r 4 )) R [Aj(-Ts),(L SK (-T 2 )Lsx(-n)) R ps]}exp(iojT 3 ) 

£ /■OO P T\ i>T 2 rT 3 

+ h dn dT2 dTz I dT4trRLsR lA'(~ T4 )' l sr{- t 3 ) Ls R {- t 2 ) L SR {- n) ps p R ] exp(iwT 4 ), 

C187 % > 
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dn dr 2 I dr 3 dr 4 (L sn L sn (- r 3 )) tr R L SR (- r 4 ) [Aj(- r 2 ) , L SK (- n) p s Pr] exp(i wr 2 ) 
Jo Jo Jo 

^ <.oo /■ti /■n /■ra 

^J ^ J ^ J J rfT4 {( isRisn ( _T4 ) isR ( _T3 ) LsR ( _ri ))R 

- (£sr£sr(-t 4 )) r (L S r(-T3)£sr(-ti)) r } [A,(-t 2 ), pg] exp(iwr 2 ) 

• roc rTi roc rT3 

+ ~R dTl / rfT2 / rfT3 / dT 4(^sRisR(-T 3 )) R (L SR (-T 4 ) J L SR (-Ti)) R [ J 4 J (-r 2 ), p s ] exp(iwT 2 ) 
ft Jo Jo Jo Jo 

^ ^CO /■Tl rT 2 PT 3 

~hj ^ J ^ J ^ J dT4 ( jLsRLsR ( _ri ))R( LsR ( _T2 ) LsR ^ T3 ))R [ A j(- T ih Ps] CXp(iWT 4 ) 
^ <.00 /■Tl 1-T2 roc 

+ kl dTl dl ~ 2 dT3 d n(LsRL SK (-T 1 )) K tr K L SR (-T 2 )[A : j(-T 3 ) , L SR (- t 4 ) p s p R ] exp(iuT 3 ) 

ft Jo JO Jo Jr 3 

^ rOO rTi rOO rT 3 

hi dTl dT2 I dl ~ 3 rfr 4 ( l sr^sr(-ti)) r [Aj(-T 2 ), (£sr(-t 4 )L sr (-t 3 )) r ps] CXp(!WT 2 ), 

(188) 

where we have performed some integral transformations as (A4). By virtue of the results of Appendixes C, E and D, 
Po^ [ = trRPiE^] (" = 2,3) can be expressed using time-integrals as 



(2) 
PQ = 



\i->+0 ' 



(3) 

Po 



(4) 

Po 



POO fT\ 

I dri / dT 2 tr R £ SR (-T 2 )L SR (-Ti)p s p R cxp(-eTi) 
Jo Jo 

rOC rT\ rT 2 

= i dTi I dr 2 / dT3tr R i SR (-T 3 )i SR (-T 2 )L SR (-Ti)p s p R exp(-eri) 
Jo Jo Jo 

rOC rT\ rT 2 rTs 

= l dTi I dr 2 / dr 3 / dr 4 tr R L SR (- r 4 ) L SK (- r 3 ) L SK (- r 2 ) X SR (- n) p s p R cxp(- e Ti) 
Jo Jo Jo Jo 



e^+0 ' 



(189) 
(190) 
(191) 



where we have considered that {TCs R (—ifi/3)) R = 0. Thus, the admittances have been expressed using the time-integrals 
alone in the fourth-order approximation for the system-reservoir interaction. In the higher order approximation, the 
admittances can be expressed using the time-integrals alone by the same way. This facilitates the calculations of 
the admittances, because when the time-correlation functions of the heat reservoir are given, the admittances can be 
calculated. I 



VI. COMPARISON WITH ANALYTICALLY 
SOLVABLE MODELS 



of the system for an external driving field which is a peri- 
odic function of the frequency u>. We adopt the following 
Hamiltonians: 



H s = huj tfb, H R = Y, huJab « b «> ( 192 ) 

a 

H SR = hbRl + hb^R, (R=^2g a b a ) (193) 



In the previous sections, we have studied expressions 
of the susceptibility in various methods. We may think 
some of them is better than others. There have been 
discussed comparative merits and demerits of the TC 
and TCL methods. Both of them are correct expression, 
and thus, if we take all the terms in expansion, they must 

be the same. When we truncate the perturbation at a where b is the boson operator of the quantum oscillator 
finite term, they are different with each other. Then, one system, b a represents the boson operator of the mode a 
of them gives the better results than the other. In this which composes the heat reservoir, uj and u a arc the 
section, we consider the two exactly solvable models, characteristic frequencies of the oscillators, and g a is the 
and compare each exact form of admittance with the coupling constant between the boson system and the bo- 
results obtained by using the methods derived in Section son of the mode a. The interaction of the quantum os- 



III in the lowest Born approximation. 
(1)^4 quantum oscillator model 



cillator system system with the external driving field is 
given by 



H ed (t) = -hbF*e iuJt ~htfFe- iut , 



(194) 



We first consider the model of a system of quantum 
oscillator interacting with a heat reservoir composed of where F is the magnitude of the external driving force. 
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(43) takes the form 

Wsr(*) = h ^2{ga bb l cx p{- i ( Lj o-u a )t} 

a 

+ g a tf 6 Q exp{i(w - uj a )i}). (195) 
The admittance x^t (w) can be written as [6] 

i f°° 

X»t(w) = ^ dtTrftfee- tLt [ft&t, p TE ]e^ 4 - e ', 

pOO ft 

= ih dtTibe- lLot cxp^{-i dTL SR (T)} 
Jo Jo 

x[b\/>rE]e iut - et , (196) 

with e— >+0. As shown in Appendix F, the above admit- 
tance can be exactly calculated as 

ih 

X66t(w) = — — (197) 

z (wo - w) + 0(w) 

with w = w + i e (e— >+0), where </>(w) is defined by 

/■OO 

#5)=/ rfr([fi(r), ijt]) R e <ffiT , (198a) 
Jo 

= *5Zlffa| 2 /(w-Wa)- (198b) 

a 

Applying the relaxation TC method (or the TCE 
method) for the Hamiltonians (192) — (194), we obtain 

tr b C (2) [w] a b f [w] = - 4>{uj) tr b o 6 t [w] + HOT, (199) 

/■ex) />r 

trbl£ ) [w} = - dr ds{(R(T)rt) K tr[hb\bp s } 
Jo Jo 

= 0, (200) 

with w = w + i e (e— >+0), where HOT denotes the higher 
order terms in powers of Hsr ■ In the lowest Born ap- 
proximation for the boson-reservoir interaction H S r, we 
have 

{i (w - w) + <j>(u>) }tr&a b t[w] = trba 6t (0), (201) 
which leads to the admittance 



x£tM = -tr hba b i[w] 



i h 



(202) 



h i (wq — w) + </>(w) ' 

with w = w + i£ (e— >+0). This is equal to the exact 
admittance (197). Therefore, for the quantum oscilla- 
tor model under consideration, the admittance (202) ob- 
tained by using the relaxation TC method (or the TCE 
method) in the lowest Born approximation for the boson- 
reservoir interaction, coincides with the exact one. 

Applying the relaxation TCL method for the Hamilto- 
nians (192) — (194), we obtain 

tr b C (2) (t) a 6t (t) = - <£(w , *) tr 6 a 6t (t) , (203) 



tr bl$\t) 



dT{(R(T)R\tv[hb\bp s ] 

-<iitfl(r)) E tr[ft6+,pB6]}e iwo(T - t) , 
= n. (9m) 



where (f>(co,t) is defined by 



0(w,t) = f dr([R(T),R^} 
Jo 



(205a) 



= i^\9 a \ 2 1 ~ ex P^^-^)^ , (205b) 



W — W Q 



In the lowest Born approximation for the boson-reservoir 
interaction Hsr, we have the TCL equation 

(d/ dt) tr b a h i (t) = —i wo tr b a b i (t) — </>(wo , t) tr b a b \ (t) , 

(206) 

which leads to the admittance [41] 



xll^) = I J o dttvnb~a b] {t)e iQt , 

= ihj dt exp ji (w — wo) t — J c£t^>(wo,t)|, 

(207) 

with w = w + i e (e^+0) . 

By using the TCLE method in the lowest Born ap- 
proximation for the boson-reservoir interaction, the ad- 
mittance was obtained by one of the authors as [41] 



TCLE / \ 
Xfefct H 



i + $(w) 



i (w - w) + 4>(coo) ' 



(208) 



with w = w + ie (e— >+0), where $(w) comes from the 
interference term and is given by 



$(w) = (<f>(w) - 0(w o ))/(w - w ). 



(209) 



Now, we have found that the admittance (202) ob- 
tained by using the relaxation TC method (or the TCE 
method) in the lowest Born approximation for the boson- 
reservoir interaction is equal to the exact one (197). On 
the other hand, the admittance (207) obtained by using 
the relaxation TCL method and the admittance (208) 
obtained by using the TCLE method take the forms dif- 
ferent from the exact one (197). 

Let us compare these admittances analytically. The 
second-order terms of these admittances in powers of the 
boson-reservoir interaction coincide and are equal to 



i he 



E 



ft I So 



(w - w ) 2 ^ (w - w ) 2 (w - ujq) 



(210) 



with u) = uj-\-i€ (e— H-0). The fourth-order terms of 



i_i _i_ 



. i_ j_ _ : l i_ 



_l* ( -i r\i-r\ F_ 
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(202)], (207) and (208) in powers of 0, as 

h<t>{oj) 2 



(4), n RTC(4), s, 



(u - (j) 3 ' 



; (w - w) 3 (uj - w a ) (w - u a >) 



TCLE(4) i \ fc 



RTCL(4) 
^b&t 



0(^0 ) 

(w - w) 3 



= V" 

^ (w - ^) 3 (W - (J a ) (W - £*><*') 



(211) 



, (212) 



E 

a , a 



h\g a \ 2 \g a '\ 2 



(a> - w) (u - uj a ) {loq - w a >) 

1 1/1 1 \ 

2 + r^TTpJ 



. (w - w) 2 

1 



+ 



2 V (w - w Q ) 2 (w — a-v) 2 

2 dJg — — W a ' 



2 ((Jo - u) ((J + (Jo - W a - (J a >) 

x(^^ + ^^)|, (213) 
\(j — (j a u — J a i / ) 

with u) = (j + i€ (e— »+0). The fourth-order term 
X^t E ' 4 '( w ) takes a form similar to the exact one x^ b \(u>). 
The fourth-order term x^^ A \(j) takes a form different 
from the exact one. It should be noted that x bb t 
and x^ L ' 4 '(o>) coincide with the exact one xj^(cj) 
in the limit lo^hjo- Thus for the quantum oscillator 
model (192) — (194), the relaxation TC method (or the 
TCE method) in the lowest Born approximation for the 
boson-reservoir interaction gives the exact admittance, 
and in that approximation the TCLE method gives the 
form similar to the exact one, while the relaxation TCL 
method gives the form different from the exact one in 
the higher-order terms. 

(2) A quantum spin model 

We next consider the model of a quantum spin system 
of magnitude S = 1/2 under an external static magnetic 
field Ho m the z direction, interacting with a heat reser- 
voir which is composed of many quantum oscillators. We 
study the susceptibility of the system for an external driv- 
ing magnetic field. We adopt the following Hamiltonians 
of the quantum spin system and heat reservoir: 



H.s = — H (Jo S z , 

H SK = -hS z (R + R t ), 



E 



h(j a blb a , (214) 
(R=Y J 9c l b a ) (215) 



where b a represents the boson operator of the mode a 
which composes the heat reservoir, u a is the character- 
istic frequency of each oscillator, and g a is the coupling 
constant between the spin system and the boson of the 
mode a. Here, uo is the Zeeman frequency uq — jH with 



quantum spin system with the external driving magnetic 
field is given by 

H ed (t) = -hS +1 ^e^ t -hS_ 1 ^e-^ t , (216) 

where H is the magnitude of the driving magnetic field. 
Then, Hsn(t) defined by the interaction representation 
(43) takes the form 

n SK (t) = -hS z Y,{9*Ke- i "° t + g* a bl e iu " «) . (217) 



In the study of this spin model, we neglect the effect 
of initial correlation between the spin and reservoir for 
brevity, and adopt the initial state p T E in the decoupled 
form 



Pte — Po Pr ; 



(218) 



with po given by (6). 

Applying the relaxation TCL method to the above spin 
model, we obtain 

tr£+C< 2 >(t)o_(t) = -V(0,t)trS+a_(t), (219) 

where 0(0, t) is defined by 

m t) = \f dr (([R(r) , i?t ] + ) R + ( [ flt( T ) , R] + ) R ), 

(220a) 



= J V|<? a | 2 (2n a + 1) /"dT(e iw - T +. 

with fi a defined by 

n a — l/{exp(/3 h(j a ) — 1}. 



— 2 UJ a T \ 



(220b) 



(221) 



Here, we have used the relations S z S + = S + /2 and 
S+S z — — S+/2 for spin S = 1/2. In the lowest Born ap- 
proximation for the spin-reservoir interaction 7Ysr, we 
have the TCL equation 



— trS+a_(i) = -ioj trS+a-(t) - ip(0,t) trS+a_(i), 



which leads the admittance 



(222) 



dttihS + a-(t) e l 



= 2ih(S z ) J dt exp |i (u — (j ) t — J dr^(0, r)|, 

(223) 

with uo = (j + i€ (e— >+0), where (S z )o is given by 
(S z )o — ^S z po- It should be noted that in the present 
spin model, the sum of the higher-order ordered cu- 
mulants vanishes, as shown in Appendix G. Therefore, 
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method in the lowest Born approximatin for the spin- 
reservoir interaction, is equal to the exact one except for 
the term of initial correlation between the spin and reser- 
voir, i.e., 



x+ _(w) = x * TCL (w). 



(224) 



Applying the relaxation TC method (or the TCE 
method) for the present spin model, we obtain 

tr S+ C {2) [w] a- [w] = - t/j(uj - w ) tr S+ 5- [w] + HOT, 

(225) 

with uj = lo + ie (e— *-+0), where ip(u>) is defined by 

1 r°° 

m=2j Q dT(([R(T),rt] + ) K + ([rt(T),R} + ) K ) 
xe Wr , (226a) 
= ij2\ 9a \ 2 (2n a + l){— ^ + ^^|. 

2 I W + W Q LU - LU a ) 

(226b) 

In the lowest Born approximation for the spin-reservoir 
interaction, we have the admittance 

x™M = itrjiS + a_M = - — 2 4^£ - v 

ft i (u>o — w) + ip{u) - W ) 

(227) 

which takes the form different from (223). 

Applying the TCLE method to the spin model under 
consideration, we obtain 

tr S+ C (2) pi [w] = - V(0) tr S +Pl [u] , (228) 

tr S+ D {2) [«] - 2 (S x ) (V(w - wo) - V(0))/(w - w ), 

(229) 

with w = w + ze (e— >+0), where tp(Q) is given by 

V>(0) = 7r^| 5Q | 2 (2n Q + l)(5K). (230) 

a 

In the lowest Born approximation for the spin-reservoir 
interaction, we have the admittance 



X T + aE (w)-2ft(S z ) 



i + ^f(u) — w ) 
i (w - w) + ip(0) ' 



(231) 



which is different from (223), where $(w — ljq) comes from 
the interference term and is given by 

*(w - w ) = (^(w - Lo ) - V(0))/(w - w ). (232) 

Thus for the quantum spin model (214) — (216), the 
admittance (223) obtained by using the relaxation TCL 
method in the lowest Born approximation for the spin- 
reservoir interaction, is equal to the exact one except 
for the term of initial correlation between the spin and 
reservoir. On the other hand, in the lowest Born ap- 



relaxation TC method (or the TCE method) and the ad- 
mittance (231) obtained by using the TCLE method have 
the forms different from the exact one. 

Let us compare these admittances analytically. The 
second-order terms of these admittances in powers of the 
spin-reservoir interaction coincide and are equal to 



2 i h (S z ) a ip(u - w )/(w - w ) 2 , 
2h(S z ) \g a \ 2 (2n a + l) 



E 



(w - W ) { (w - W ) 2 - w 2 } ' 



(233) 



The fourth-order terms of those admittances are ob- 
tained, by expanding (223), (227) and (231) in powers 
of tp, as 



/■oo />£ />t 

X+ C - {i) (u)=th(S z } / dt / dn / dr, 
Jo Jo Jo 



(ir> - 



x V(0,Ti)V>(0,T 2 )exp{i(w-w )i}, 
ft (Sz)o I g a I 2 • I g a > | 2 (2 n a + 1) (2 ri Q < + 1) 



J, (w - w ){(w - w ) 2 - w 2 }{(w - w ) 2 - w 2 ,} 

x 6 (w - wp) 4 - 3 (w - wp) 2 (w 2 + w 2 ,) + (w 2 - w 2 ,) 2 
{(w - w ) 2 - (w Q + w Q /) 2 }{(w - w ) 2 - (w Q - W Q /) 2 } ' 

(234) 

X * TC(4) (w) - 2 ft (S z ) V(w - w ) 2 /(w - w ) 3 , 

= - 2 ft (S z ) | g a | 2 ■ | g a > | 2 (2 n a + 1)(2 fW + 1) 
(w - ^o){(w - w ) 2 - w 2 }{(w - Wo )2 - W 2 ,} ' 

(235) 

x TCLE(4) (w) = 2 fi ( ^ )o _ ?/)(5)/( _ _ ^ )3j 

2 i Tr ft (&) | g a | 2 ■ | g a / | 2 (2 n a + 1)(2 fW + 1) 
(w-w ) 2 {(w-w ) 2 -w 2 } 

xj(w a ,), (236) 



E 



with w = w + ze (e— >+0). These fourth-order terms of 
admittances take the forms different from each other. 
They do not coincide even in the limit w^w . For this 
quantum spin model, the relaxation TC method (or the 
TCE method) and the TCLE method in the lowest Born 
approximation for the spin-reservoir interaction give the 
forms of admittance different from the exact one in the 
higher-order terms. 

For the quantum oscillator model (192) — (194), the re- 
laxation TC method (or the TCE method) in the low- 
est Born approximation for the boson-reservoir interac- 
tion gives the exact admittance, and in that approxima- 
tion the relaxation TCL method and the TCLE method 
give the forms of admittance different from the exact one 
in the higher-order terms. On the other hand, for the 
quantum spin model (214) — (216), the relaxation TCL 
method in the lowest Born approximation for the spin- 
reservoir interaction gives the exact admittance except 
for the term of initial correlation between the spin and 
reservoir, and in that approximation the relaxation TC 
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the forms of admittance different from the exact one in 
the higher-order terms. This means that accuracy of the 
admittances obtained by using the relaxation TC method 
(or the TCE method), the relaxation TCL method and 
the TCLE method in the lowest Born approximation for 
the system-reservoir interaction, depends on the model 
of physical system under consideration. 



VII. SUMMARY AND CONCLUDING 
REMARKS 

We have derived the exact expressions (38), (56), (72) 
and (91) of the complex admittance for quantum sys- 
tem in contact with heat reservoir using the relaxation 
method and the external-field method, respectively, for 
the two types of equation of motion, which are the 
time-convolution (TC) equation and time-convolutionlcss 
(TCL) equation. We have shown that the expression 
of the complex admittance obtained using the relax- 
ation method with the TC equation (the relaxation TC 
method), coincides exactly with that obtained using the 
external-field method with the TC equation (the TCE 
method), i.e., (99). We have also shown that the three 
expressions of the complex admittance coincide with each 
other in the lowest Born approximation for the system- 
reservoir interaction, i.e., (127), though the three expres- 
sions of the exact admittances are different from each 
other. We have besides given the expressions of the ad- 
mittance in the n-th order approximation for the system- 
reservoir interaction, and have expressed the admittance 
in the fourth order approximation using time-integrals 
alone by transforming inverse-temperature-integrals into 
time-integrals. Furthermore, we have compared the ad- 
mittances for the two exactly solvable models analyti- 
cally and have found that for the quantum oscillator 
model (192)-(194), the relaxation TC method (or the 
TCE method) in the lowest Born approximation for the 
boson-reservoir interaction gives the exact admittance, 
that for the quantum spin model (214) — (216), the relax- 
ation TCL method in the lowest Born approximation for 
the spin-reservoir interaction gives the exact admittance 
except for the term of initial correlation between the spin 
and reservoir, and thus that accuracy of the admittances 
depends on the model of physical system under consid- 
eration. 

We here compare the admittances (44) [or (77)], (60) 
and (97), which are, respectively, obtained using the re- 
laxation TC method [or the TCE method], the relax- 
ation TCL method and the TCLE method in the lowest 
Born approximation for the system-reservoir interaction, 
with the results in the conventional Markovian approx- 
imation [27, 29-32] or in the van Hove limit [51]. In 
the conventional Markovian approximation, the motion 
of the system is determined in the van Hove limit [51] or 
in the narrowing limit [27, 29] in which the heat reser- 
voir is damped very rapidly, i.e., the reservoir correlation 



(40) and the second-order TCL equation (57) become 
[27, 29, 30, 41, 43] 

(d/dt) a(t) = -iL s a(t) + C {2) a(t) , (237) 

where is given by (94). Substituting the formal so- 
lution of the above equation : 

a(t) = cxp{- i L s t + C (2) t}d,(0) (238) 

into the Kubo formula (16), the admittance takes the 
form [41-43] 

Xj» - (i/h) tvA^i (L s -w)- C^y 1 [A 3 , Ps + p^] , 

(239) 

which is valid only in the van Hove limit [51] or in the 
narrowing limit [27, 29-32], where p is given by (45) or 
(189). Comparing the above admittance (239) with the 
admittance (97) obtained using the TCLE method in the 
lowest Born approximation for the system-reservoir in- 
teraction, the latter admittance includes the interference 
term Dj [cj] which is not included in the former admit- 
tance and represents the effects of the initial correlation 
and memory for the system and reservoir. The effects 
of the initial correlation and memory are the effects of 
the deviation from the van Hove limit [51] or the nar- 
rowing limit, because the effects of initial correlation can 
be neglected in that limit as shown in Ref. [42] and the 
memory effects are the effects of collision of the system 
with the heat reservoir. These effects are neglected in 
the conventional Markovian approximation [27, 29-32]. 
As examined in Section IV, the admittance obtained us- 
ing the relaxation TC method [or the TCE method], the 
admittance obtained using the relaxation TCL method 
and the admittance obtained using the TCLE method, 
coincide with each other in the lowest Born approxima- 
tion for the system-reservoir interaction. Therefore, the 
admittances (44) [or (77)], (60) and (97) include the ef- 
fects of the initial correlation and memory for the system 
and reservoir. 

The relaxation method, in which the Kubo formula is 
calculated for systems with no external driving fields, is 
essentially equal with the external-field method in which 
the admittance is directly calculated from equations of 
motion with external driving terms. Because, when the 
TC equation of motion is used to calculate the admit- 
tance, the expression of the complex admittance obtained 
using the relaxation TC method, coincides exactly with 
that obtained using the TCE method (the external-field 
method with the TC equation). When the TCL equa- 
tion of motion is used to calculate the admittance, the 
two expressions of the complex admittance obtained us- 
ing the relaxation TCL method and TCLE method (the 
external-field method with the TCL equation) arc differ- 
ent to each other, and also they arc different from the 
expression of the admittance obtained using the relax- 
ation TC method (or the TCE method). The difference 
of the expressions occurs by renormalizing the memory 
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In Section IV, we have examined analytically the rela- 
tions among the admittances (44) [or (77)], (60) and (97), 
which are obtained using the relaxation TC method [or 
the TCE method], the relaxation TCL method and the 
TCLE method, respectively, in the lowest Born approx- 
imation for the system-reservoir interaction. The fact 
that the three expressions (44) [or (77)], (60) and (97) 
of the complex admittance coincide with each other in 
the lowest Born approximation for the system-reservoir 
interaction, means that the three expressions of the ad- 
mittance are respectively valid in that approximation. 
The differences among the three expressions come from 
the fact that the high-order terms differ from each other. 
These three expressions include the high-order terms that 
are dominant in the resonance region. Even if the expres- 
sion of the complex admittance obtained in the lowest 
Born approximation for the system-reservoir interaction 
includes many more high-order terms, it is not always the 
more exact expression. Because, that expression of the 
complex admittance does not include part of the third- 
order, fourth-order and high-order terms in powers of 
the system-reservoir interaction. It is also possible to de- 
rive the admittance (60) obtained using the relaxation 
TCL method from the admittance (97) obtained using 
the TCLE method by proceeding in the same way as in 
Section IV. Moreover, it should be noticed that as men- 
tioned in Section III, the admittance (44) obtained using 
the relaxation TC method by the lowest Born approxi- 
mation of the equation (28) including no external driv- 
ing terms, coincides with the admittance (77) obtained 
using the TCE method by the lowest Born approxima- 
tion of the equation (63) or (68) including external driv- 
ing terms. Thus in the relaxation TC method and TCE 
method, although the equations in which the perturba- 
tion is truncated are different from each other, the ex- 
pressions of the obtained admittances coincides with each 
other. 

The three expressions (44) [or (77)], (60) and (97) of 
the complex admittance are respectively valid in the low- 
est Born approximation for the system-reservoir inter- 
action as mentioned in the above paragraph, but the 
usability is different from each other. The expression 
(60), which is obtained using the relaxation TCL method, 
takes the form that is not easy to use, because it is dif- 
ficult in general to calculate the time integral of the or- 
dered exponential for many-body systems or for com- 
plicated systems. The expression (97), which is ob- 
tained using the TCLE method, has been formulated in 
terms of thermo-field dynamics (TFD) using the method 
of noncquilibrium thermo-field dynamics proposed by 
Arimitsu and Umezawa [52] in order to study the lin- 
ear response of many-body systems interacting with heat 
reservoir [43, 44, 47]. It has been applied to an weakly- 
interacting boson system [53] and to a ferromagnetic sys- 
tem interacting with a phonon reservoir in the spin- wave 
approximation [54] , by using the generalized thermo-field 
dynamics [43, 44, 47]. Recently, fluctuation-dissipation 
theorem for the TCLE method has been examined in the 



lowest Born approximation for the system-reservoir inter- 
action [47, 53]. The TCLE method is a method for im- 
proving the difficulty of the relaxation TCL method. The 
expression (44) [or (77)], which is obtained using the re- 
laxation TC method [or the TCE method] , can be applied 
to a interacting spin system by making a unfied numeri- 
cal method [37]. In practical applications, the relaxation 
TC method, the TCE method and the TCLE method 
may be the useful methods to study the linear response 
of many-body systems interacting with heat reservoir. 

We have derived the formulae necessary for the higher 
order expansions in powers of the system-reservoir in- 
teraction in Section V. When we proceed to the high 
order in powers of the system-reservoir interaction, the 
formulae derived in Section V will be useful. We have 
there given the expressions of the admittance in the n- 
th order approximation, and have also given the forms 
of the admittances in the fourth-order approximation in 
the expression using time-integrals alone by transforming 
inverse-temperature-integrals into time-integrals. This 
facilitates the calculations of the admittances, because 
when the time-correlation functions of the heat reservoir 
are given, the admittances can be calculated. 

As examined in Section VI for the two exactly solvable 
models analytically, the relaxation TC method (or the 
TCE method) in the lowest Born approximation in pow- 
ers of the system-reservoir interaction gives the exact ad- 
mittance for the quantum oscillator model (192) — (194), 
and the relaxation TCL method in the lowest Born ap- 
proximation in powers of the system-reservoir interac- 
tion gives the exact admittance for the quantum spin 
model (214) — (216) except for the term of initial corre- 
lation between the spin and reservoir. This means that 
accuracy of the admittances obtained using the relax- 
ation TC method (or the TCE method), the relaxation 
TCL method and the TCLE method in the lowest Born 
approximation for the system-reservoir interaction, de- 
pends on the model of physical system under considera- 
tion. For the case of a non-adiabatic interaction as (193), 
in which the Hamiltonian 7i s of the quantum system is 
non-commutable to the system-reservoir interaction H.sr, 
the relaxation TC method (or the TCE method) may give 
a better result. For the case of an adiabatic interaction 
as (215), in which the Hamiltonian TL S of the quantum 
system is commutable to the system-reservoir interaction 
Hsr, the relaxation TCL method may give a better re- 
sult. Accuracy of the three admittances is an interesting 
subject of future study. 
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APPENDIX A: TRANSFORMATION OF THE /^'-INTEGRAL INTO A TIME-INTEGRAL 

We show that the expression (74) for the interference term [u>] in the TC equation, which is equal to the first term 
of the interference term Z)( 2 )[w] given by (95), coincides with the expression (75). The inverse-temperature-integral 
(/^'-integral) in the expression (74) can be integrated and can be rewritten as follows, 

/ d(3'cxp(-iL Q T)L ed [uj}p s p R H SK (-ih(3')e^ T = / dp cxp(- i L r)L ed [iu} p s p K exp(/3' H L ) H SR e l " T , (Al) 
Jo Jo 

= cxp ( - i L Q t) L ed [u] p s p R cx P(^ fi Z/ °) L WgR e l0JT =exp(-iL r) L ed [w] ^- L SK p s p R e l u T , ( A2) 



% U) T — £ Z 



(A3) 



which can be expressed using a time-integrals as 

/ d(3' e- lLoT L ed [uj}p s p K H S R(~ihp')e luJT = lim i ds e~ 1 L ° T L ed [uj] e l L ° (r " s) L SR p s p R e 
Jo e ^+° Jt 

Using the integral transformation 

/>oo />oo poo ps roc rT 

/ dr ds= ds dr = dr ds(r^s), (A4) 

Jo Jt Jo Jo Jo Jo 

the expression (74) for ZM 2 )[o>] coincides with the expression (75). Similarly, the expression (47) for D^[uj}, which is 

(2) 

equal to the first term of D v - '[uj] given by (98), can be shown to coincide with the expression (78). 

APPENDIX B: DERIVATION OF THE INTERFERENCE TERM D[u>] 
The interference term D(t) given by (83), in the limit t — > oo, becomes 



D(t)=[ dTti K L SK 1 e- lQLQT \-L ed (t-T)Qp TE + Z(t-r) L ed (t - r) p TE } I 
Jo ^ 1 - M J J 



(Bl) 



(B2) 



where we have substituted (8) and E(t — r), which is given by 

Z(t-T) = -ij dscxp{-iQLQ(s-T)}QL SR Taq>{iL{s-T)}. (B3) 
Then, the interference term D[cu] takes, according to the definition (87), the form 

,oo 1 

D[uj]=- I drtr R L S R-; ^- exp(- i Q L Qt) L ed {uj] Q p TE cxp(iu; t) 

Jo 1 — 2j 

/•OO f-OO -j^ 

-i J dr J dstr R L SR y— — cxp(- i Q L Qs) QL SK P cxp{i L (s - T)}L ed [Lu] p TE exp(iuj t), (B4) 
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APPENDIX C: TRANSFORMATION OF THE ft AND ft INTEGRALS INTO TWO TIME-INTEGRALS 



The ft and ft integrals in Subsection V.C can be integrated and can be rewritten as follows, 
rP fPi 

/ dft / dfoe-* L ° T [A^spzHsni-ihp^Hsni-ihfo)], 
Jo Jo 

= f d$ x f 1 d(3 2 e- iL °T[A,p sPK (e^ hL °H sn )(e^ hL °H sn )], 
Jo Jo 

= f dfr e~ 1 LaT \A , p s p R { h L ° H SR -±r- H SR - (e^ 1 h L « H SR ) ^— H sl 
Jo 1 hL 'hL a 

f ePKLo-1 s i 



-iLnT 



= e- lLoT 



e PhL _ -y i , e f)hL () _ -y 

A ' hT [ WsR ^ WsR ' Ps PR ] ~ (hT [ HsR ' Ps PR ] ) ^ HsR 



(CI) 
(C2) 
(C3) 
(C4) 
(C5) 



which can be expressed using the two time-integrals as 
f P f Pi 

/ dft / d/3 2 e- 4ioT [A, PsPnHszi-ih^Hsni-ihfo)], 
Jo Jo 



i 



dri J dT 2 e- lLoT [A, e 



i L Q (n - t) 



[w 



SRC 



- i L a (t 2 - Ti) 



H S R , PS PR ] _ 



IT1 
POO 



I £-►+() 



. (C6) 



dn / dr 2 e- 4 io T [A , (e" 1 Lo ^ " ^ [ H SK , p s p R } ) e~ 1 L « ^ ~ T > H SK ] } exp(- e n - e r 2 ) | 
= ^ { y dn y dr 2 [A(- r) , [ W SR (- n) W SR (- r 2 ) , p s p R ] ] 

/OO /"GO 
dn J dT 2 [A(-T),[H SR (-n),p s p K ]H SR (-T 2 )]}exp(-eT 1 -eT 2 )\^ +0 , (C7) 

= ^2 y dTl y rfT2 l/ 4 ( _T ) ' ^sr(-T2)[W sr (-ti) , p s p R ] - [Hsr(-ti) , p s Pr ] H SR (- t 2 ) ] ] exp(-£T 1 )| £ ^ +0 , 

■2 roc rn 

= jpj d Tl J dr 2 [A(- T ),[H SR (-T 2 ), [W S R(-r 1 ),p s p R ]]]exp(- e r 1 )|^ +0 , (C8) 
dn J dr 2 [A(-t), L sr (-t 2 )Lsr(-ti)PsPr] exp(-en)| e ^ +0 - (C9) 



APPENDIX D: TRANSFORMATION OF THE ft, ft AND ft INTEGRALS INTO THREE 

TIME- INTEGRALS 

The ft, ft and ft integrals in Subsection V.C can be integrated and can be rewritten as follows, 

f dft f 1 dp 2 l h dft e~ lL ° T [A, p s Pr Wsr(- ifiA) W SR (- i fift) W S r(- ift/?s) ] , 
Jo Jo Jo 

e /3i RL _ ^ 



d(3i e 



-iLnT 




1 

fi3 



fi3 



, „/3i RLo _ 1 1 

^4 , p s Pr (e 01 fi io W SR ) { r7 W SR — H SR 



1 



i L a t 



i La t 



7 ^0 -^0 



(Dl) 
(D2) 



-4, PsPR-fy- (e fihLa - 1) Hsr ^ Hsr ^ W S r - (-!- (e 

^0 ^0 -^0 v -^0 

- f T- ( e/3RL ° - 1) ^sr ^ Hsr) ^ W S r + f-i- {e fihLa - l) W SR ) (-1- H S r) ^- H S r}} , (D3) 

-4 , — [ W S r — Hsr ^— Wsr , Ps Pr — f ^— [ Wsr , Ps Pr ] \ ~r~ Wsr -7— Wsr 

_Lo L -^0 -^0 J Vi 7 -^0 L o 

- [n SR -^H S R, PsPr]) ^H S r+ [W S r, PsPr])(^W S r) ^HsrI, (D4) 
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which can be expressed using the three time- integrals, by performing some integral transformations as (A4), (138) 
and (142), as 

rP 



f dfr [ 1 d/3 2 [ 2 d/3 3 e- iLaT [a, Ps p K n SK (-ihp 1 )n S n(-ihp 2 )n S n(-ihp 3 )], 

Jo Jo Jo 



^3 r foo poo 



Tl JT 2 

oc /■ :x: 



r Jt 2 

roc poo 



[•OO [OO 



dr 3 e- iLoT [A, e~ 


i L (n - 




- i Lo (r 2 - .0 e - i L (r 3 - r 2 ) ^ , pg pR ] ] 


dr 3 e-' LaT [A, (e 


-iL (n 


- T) [«SR, 


Ps Pr ] ) e" 4 L ° (T2 - T) W SR e" 4 io (T3 - T2) Hsk ] 


dr 3 e- lLoT [A, (e 


- i L (ti 


~ r) [Hsae 


-,Lo ^-^) WsR)/0si0R ]) e -iic (-3— ) WsR ] 


dT 3 e~* L ° T [A, (e" 


- i L (n 




PsPr]) {e- lLo (T2 - T) W S R)e- 4io(T3 - r) H S R]} 



X CX P (-£T 1 -£T 2 -£T 3 )| £ ^ +0 , (D5) 

j3 , /-oo |-oo /-oo 

= ¥lj dTl J dT2 J dT3 r) ' [ Hsr(_ Tl) Wsr( ~ T2) Hsr( ~ T3) "° s Pr ] ] 

/■ OO p oo ^ oo 

- / dn dr 2 / dr 3 [A(- t) , [ W SR (- n) , p s Pr] Wse(- Ti) Wsh(- t 3 ) ] 

/>00 />OC />OC 

-/ dn / dr 2 / dr 3 [A(-r) , [H sr {-ti) H S r(- t 2 ) , PsPr] H S r(-t 3 )] 

J T J Tl J T 

/OO />OC /■ OO 

dn y dr 2 y dr 3 [^(-^^^(-n), PsPR]H S R(-r 2 )H SR (-r 3 )]}e- eTl - £T2 - eT3 |^ +0 , (D6) 

•3 /-OO /.Tl />T2 

= ^y dn y dr 2 y dr 3 [A(-T),[W S R(-T3),[W S R(-T 2 ),[W S R(-n),PsPR]]]]exp(-en)|^ +0 , (D7) 

/OO /*Ti pT2 

dn J dr 2 y dr 3 [A(-r), J L SR (-r 3 )L SR (-r 2 ) J L SR (-Ti)p s/ 9R] exp(-en)| e ^ +0 - (D8) 



APPENDIX E: TRANSFORMATION OF THE fa, fo, fa AND /? 4 INTEGRALS INTO FOUR 

TIME- INTEGRALS 

The /3i, /3 2 , /3 3 and /3 4 integrals in the expression (159) of p$ can be integrated and can be expressed using the 
four time-integrals by proceeding in the same way as in Appendixes C and D as follows, 

rP rPi rP? rP* 



( dp! f 1 d(3 2 f 2 d(3 3 f 3 d/3 4 psPRH SR (-2ft/3i)H SR (-2ft/3 2 )H SR (-2ft/3 3 )H SR (-ift/3 4 ), (El) 
Jo Jo Jo Jo 

± f d/?i Ps p R ( e ft h La Hsk ) f f (e^ h L ° l) H S r f H S k f W SR - (f (eft - l) W SR ) f W SR f W SR 
ft Jo ^0 -^0 ^0 v -^0 7 ^0 L o 

(±_ ( e ft _ 1} HsR J. Wg \ J. HsR + (■ J. ( e ft HL 1} HgR \ / I n \ 1 Wsr \ (E2) 

i Ps PR {y- {e' 3hLa - 1) Hsr f H SR f H SR f H SR - (f (e f3hL ° - l) W SR ) f H SK f W SR f H SR 
- f T- ( e/3 " i0 - 1) ^sr f W SR ) f Hsr y- Hsk + (-- {e f3hL ° - l) Hsk) (-- Hsk) f Hsk t~ Hsk 

\Lo Lq / Lo -Lq \Lo / \Lq / Lq Lq 

_ (_L ( e /J"o _ 1) HsR J_ WsR J_ U \ J_ WgR + C J_ ( e /3^c _ !) HsR ^) (J_ HsR J. HsR ^) J_ WsR 

+ (J- ( e ^io _ i) HsR J_ WsR ) (_L HsR ) J_ WsR _ ( J_( e /3^o _ !) WgR ) (_L HsR ) C J_ WgR ) J. WsR \ (E3) 
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i 

¥ 



dri / dr 2 / dr 3 / dr 4 [ W SR (- n) H SK (- r 2 ) H S r(- r 3 ) W SR (- ^4) , Ps Pr ] 



r 1 >^ T2 ^ T3 

00 />oo 



/>oo />oo />oo />oo 

/ dri I dr 2 I dr 3 / dr 4 

JO JO Jt 2 Jt 3 



I 

¥ 



[Wsr(-ti) , PsPr] H S r(- t 2 ) H S r(- t 3 ) H S r(- t 4 ) 
[Wsr(-ti) W sr (-t 2 ) , p s p R ] H S r(-t 3 )Hsr(-t 4 ) 
[Wsr(-ti) , PsPr] H S r(-T2)Hsr(-T3)Hsr(-t 4 ) 
[Wsr(-ti) W sr (-t 2 ) W S r(-t 3 ) , PsPr] W S r(-t 4 ) 
[Wsr(-ti) , PsPr] H S r(-T2)Hsr(-t 3 )Hsr(-t 4 ) 
[Wsr(-ti) W SR (-r 2 ) , PsPr] Hsr(-t 3 )Hsr(-t 4 ) 

[Wsr(-Ti), p S p R ] HsR(-T 2 )HsR(-T 3 )?i:sR(-T 4 )|e _ 

dn I dr 2 I dr 3 I dr 4 [ H S r(- n ) H S r(- t 2 ) H S r(- r 3 ) W SR (- 7" 4 ) , Ps Pr ] 

/•OO /*CO /*oo /*CO 

- / dn / dr 2 / dr 3 / dr 4 [ W SR (- n) W SR (- r 2 ) W SR (- r 3 ) , p s Pr ] W SR (- r 4 ) 

/•OO />00 />00 /-T3 

+ / dri / dr 2 / dr 3 / dr 4 [ W SR (- n) W SR (- t 2 ) , p s Pr] W SR (- r 3 ) W S r(- t 4 ) 
Jo Jti jo jo 

/•OO />00 />T2 ^3 

/ dri/ dr 2 / dr 3 / dr 4 [ W SR (- n) , p s p R ] W SR (- r 2 ) W SR (- r 3 ) W SR (- r 4 ) | 
Jo Jo Jo Jo ' 



1-3 

/* 00 /* 00 />oo 

dn / dr 2 / dr 3 / dr 4 

Jti JO Jt 3 

/>00 ^OO ^OO /'OO 

+ / dri / dr 2 / dr 3 / dr 4 

JO JO JO Jt 3 

/•oo /* 00 /* 00 />oo 

— / dri / dr 2 / dr 3 / dr 4 

JO Jti Jt 2 Jo 

/>oo />00 />00 />00 

+ / dri / dr 2 / dr 3 / dr 4 

JO JO Jr 2 JO 

/•OO i° OO /*00 /'OO 

+ / dn / dr 2 / dr 3 / dr 4 

Jo Jti Jo Jo 

/•OO /-OO />00 /* OO 

/ dri / dr 2 / dr 3 / dr 4 

Jo Jo Jo Jo 



e Ti— e T2 — e t 3 — e t 4 



le^+0 ' 

(E4) 



x exp{- eri - er 2 - er 3 - er 4 }| e ^ 



e^+o ' 



(E5) 



which can be expressed, by performing some integral transformations as (A4), (138) and (142), in the following 
compact form 

f * dfc f 1 dfo f 2 d(3 3 / /33 d/3 4 p s p R 7^ SR (-^ft/3 1 )7^ SR (-^^/3 2 )H SR (-^fi/3 3 )WsR(-^/34), 
Jo Jo Jo Jo 

^4 P OO I-Ti /.T 2 /-T 3 

= ¥j dTl J dT2 J Q dT3 J Q dr 4[ W sR(-r 4 ),[7i SR (-r 3 ),[W SR (-r 2 ),[W SR (-r 1 ),p s p R ]]]]e- eTl |^+o, (E6) 

f' OO /'Ti />T2 ^3 

= z 4 / dn / dr 2 / dr 3 / dr 4 L S r(- r 4 ) X SR (- r 3 ) L S r(- r 2 ) L S r(- n) Ps Pr exp(- e n) | e ^ +0 . (E7) 
Jo Jo Jo Jo 



APPENDIX F: CALCULATION OF ADMITTANCE Xw,t(w) [(196)] 

The admittance X6fet(w) given by (196) can be rewritten as 

/•oo r 00 f Tl ri~ n -i ^ 

X w(uj)=ifr dtTr6(t){l + V(-z) n / dn / dr 2 • • • / dr„ L SR (n) L SR (r 2 ) • • • i SR (r„) [ 
Jo L n=1 Jo Jo Jo ' 

x [b\ p TE ]cxp(iujt-et) | e ^ +Q , (Fl) 
= ih\ dtTr6+ V(-z)" / dr„ / dr„_i • • • / dn/ dt Tr [•••[[ & , W SR (n) ], H SR (r 2 ) ],•• • 

L JO n=1 Jo Jr n Jt 2 Jti 

••• , n S R(T n )]/h n ) \b\ pj E }expU(ij-u; )t-et}\^ n , (F2) 
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which can be calculated using (195) as 

X66t(w) = rH rVl-sf" 1 / dr 2m dT 2m -i ■ ■ ■ / dr 2 / rfn V • • • V | g Ql | 2 • • • | g, 



2 



m=l " U '2m "T 3 ^T 2 Ql Qm 

x Tr b [ 6 f , /9 TE ] exp{i (u> - u ) t\ + i (u> - u> ai ) n - i (lj q - uj ai ) t 2 H 

hi (w - w Qm )r 2m -i - i (w - w am )T 2m }, (F3) 

= ^( 1+ fy v (-ir\g ai \ 2 ---\g am \ 2 (m 

wo-wl £^ (wo-w) m (w-w Ql ) ••• (w-w ara ) ' 



fi + y l^d r= if , ( F5) 



wo ^ w ^ (wo — w) (w — w ) J i (wo — w) + </>(w) 

with w = w + ie (e— >+0), where </>(a>) is given by (198). Thus, the exact form (197) of Xbwi^) can be obtained. 



APPENDIX G: THE ORDERED CUMULANTS FOR THE SPIN MODEL (214) (216) 

We calculate the ordered cumulants for the quantum spin model (214) — (216). For the second-order ordered cumu- 
lant {L SR L SK (— t)) r (t > 0), we can obtain, by using the relations S Z S+ = S+/2 and S+S z = — S+/2, 

tvS + (L sn L SR (-T)) K a_(t) = (l/2){([R(r), Jjt] + ) a + ( [i?t( r ), R] + ) R } tr S+ o_(t). (Gl) 

We can calculate the fourth-order ordered cumulant (L SR L SR (— t 3 )L sr (— t 2 )L sr (— t i))r ( t i > t" 2 > t 3 > 0) as follows, 

tr S+ (L SR L SR (- r 3 ) L SR (- r 2 ) L SR (- n)) R 5_ (t) = - tr [S+, S z } ( (R + i? f ) L SR (- 73) Lsb(- T a) £sr(- n) )r «_(*)> 
= - Tr {5+ S z (i2 + i? f )(i?(- r 3 ) + J2t(- 73)) - S z S+ (R(- r 3 ) + flt(- t 3 ))(# + fit)} L SR (- r 2 ) L SR (- n) ^a_(t), 
= (1/2) Tr 5+ [ i? + fit, T3 ) + fl t(_ T3 ) ] + LsR (_ T2 ) isR (_ Tl ) pR a_ ( t ) f 

= (1/ 2 2 ) Tr S+ [ [ J2 + R\ R{- 73) + i? f (- r 3 ) ]+ , r 2 ) + flt(- r 2 ) ]+ L SR (- n) p R 5_ (t), 

= (1/ 2 3 ) ( [ [ [R + R\ R(- r 3 ) + R\- r 3 ) ]+ , R(- r 2 ) + flt(- r 2 ) ]+ , n) + flt(- n ) ]+ ) R tr 5+ 5_(t), (G2) 

which can be rewritten using the second-order ordered cumulants by the Wick's theorem for finite temperature, as 

tr S+ (Lsk L S n{- r 3 ) L SR (- t 2 ) L sr (- n)) R a_(i), 
= (1/ 2 2 ) {( [R + R\ R(- r 3 ) + i?t(- T3 ) ]+), ( [R(- r 2 ) + flt(- r 2 ) , R(- n ) + i? f (- n) ]+) R 
+ ( [fl + i2t, T 2 ) + flt(- T 2 ) ] + ) R ( r 3 ) + i?t(- r 3 ) , Tl) + iZ f (- Tl) ]+> R 

+ ([R + R\ R(- ti) + i?t(- ti) ] + ) R ( [ J2(- t 3 ) + i? f (- r 3 ) , t 2 ) + J2t(- t 2 ) ] + ) R } tr S+ 5_ (t). (G3) 
We can also calculate the product of the second-order ordered cumulants as 
tr S+ (L S r ^sr(- t 3 )) r (L S r(- t 2 ) L sr (- ti)) r a_ (i), 

= (l/2)([i2 + i2t, iZ(-T 3 ) + iZt(-T 3 )] + ) R trS+(L SR (-T 2 )L SR (-Ti)) R a_(t), 

= (1/ 2 2 ) ( [R + i?t, R(- T3 ) + i?t(- T3 ) ]+) E ( t 2 ) + RH- t 2 ) , ti) + flt(- ti) ]+) E tr 5+ 5_(t), (G4) 
tr S* + (L SR L S r(- t 2 ))r (i S R(- 7-3) L SR (- ti)) r a_ (i), 

= (1/ 2 2 ) ( [i? + i?t, i?(- T2 ) + i?t(- r 2 ) ]+) B ( t 3 ) + flt(- t 3 ) , i?(- ti) + i?t(- ti) ]+) B tr 5+ 5_(t), (G5) 
tr S* + (L SK L SR (- ti)) r (L S r(- 73) £sb(- t 2 ))r a_ (i), 

= (1/ 2 2 ) ( [i? + i?t, J2(- ti) + i?t(- ti) ]+) B ([R(- t 3 ) + i? f (- r 3 ) , i?(- t 2 ) + i?t(- t 2 ) ]+) B tr 5+ 5_(i). (G6) 

Therefore, the sum of the fourth-order ordered cumulants vanishes, i.e., 

tr S* + {(L SR L sr (- t 3 )L sr (- t 2 )L sr (- ti)) r - (L SR L sr (- r 3 )) R (L SR (- t 2 )L sr (- ti)) r 

- (£sb Lsr(- t 2 )) r (L S r(- t 3 )L sr (- ti)) r - (L SR L SR (- ri)) R ( J L SR (- t 3 )L sr (- t 2 )) r } 5_ (t) - 0. (G7) 



r - i_ • _i- 
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